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1.  INTRODUCTION 


The  problem  of  recognizing  patterns  in  a  two-dimensional 
image  is  an  extremely  general  one,  and  we  must  break  it  up  into 
more  limited  problems  in  order  to  make  progress.  We  know  from 
the  transmission  of  information  that  we  can  divide  the  problem 
of  reception  of  time  signals  -  or  pattern  recognition  with  one 
variable  -  into  two  parts.  First,  we  have  the  reception  of 
one  signal  out  of  a  small  set  of  possible  signals.  This  problem 
is  encountered,  e.g.,  in  the  reception  of  teletype  signals.  We 
know,  due  to  the  synchronization  between  transmitter  and  receiver, 
when  a  signal  is  received,  and  we  also  know  that  it  must  be  one 
of  the  32  possible  teletype  signals.  The  preferred  way  to  decide 
which  signal  was  actually  transmitted  is  to  multiply  the  received 
signal  with  samples  of  the  32  possible  signals,  and  to  integrate 
the  32  products.  The  sample  signal  producing  the  largest  value 
of  the  integral  is  accepted  to  be  the  one  which  equals  the  trans¬ 
mitted  signal.  This  process  is  usually  referred  to  as  signal 
detection  by  means  of  cross-correlation  with  sample  signals. 

Very  different  is  the  detection  of  the  (time)  patterns  of 
telephone  signals.  Even  though  the  number  of  telephone  signals 
transmittable  during  any  given  time  interval  is  finite,  the 
number  is  so  large  that  cross-correlation  with  sample  functions 
becomes  impractical.  The  ear  uses  a  detection  process  that  may 
be  described  as  a  sliding  correlation  based  on  the  Fourier 
series . 

We  can  apply  the  same  classification  to  two-dimensional 
(spatial)  patterns  with  the  variables  x  and  y  rather  than  the 
one  time  variable  t.  A  machine  that  can  detect  the  characters 
of  a  typewriter  -  jf  which  there  are  about  100  -  needs  to 

be  "synchroni zed"  to  the  location  of  the  characters,  and  may  then 
use  cross-correlation  with  100  sample  characters.  This  process 
cannot  be  used  to  detect  handwritten  characters,  or  patterns  on 
aerial  photographs,  since  the  number  of  possible  "sample  charac¬ 
ters"  would  be  enormous,  even  though  it  would  be  finite.  We 
need  something  comparable  to  the  detection  of  audio  signals  by  the 
ears . 
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Lot  us  look  at  I-'ig.  1.  We  recognize  instantly  several 
triangles  at  various  locations  with  various  rotations,  sizes, 
and  representations,  hut  equal  shape.  With  representations 
we  mean  that  a  triangle  may  be  represented  by  a  black  line  with 
three  corners,  by  a  black  area  with  three  edges,  and  in  many 
other  ways.  Hqual  shape  means  that  the  triangles  can  be  made 
to  coincide  by  shifting,  rotating,  changing  the  size,  and  chang¬ 
ing  the  representation.  Alternately,  two  triangles  with 
angles  ctj  ,  Vj  and  c. , ,  y,  -  listed  in  both  cases  either 

clockwise  or  anti  clockwise  -  have  the  same  shape  for  aj=cs, 

Pl"P:.  vi'V 

If  we  had  used  correlation  with  sample  triangles  to  decide 
that  I-ig.  1  shows  triangles,  we  would  have  to  correlate  triangle 
with  various  shapes,  representations,  x-shifts,  y-shifts,  rota¬ 
tions  and  sizes.  This  would  require  so  much  time  that  we  still 
would  not  have  come  to  a  decision.  Obviously,  our  brain  does 
not  use  correlation  to  recognize  a  triangle. 

How  do  we  recognize  a  triangle?  The  word  implies  that  the 
three  angles  are  the  important  feature.  The  angles  are  indepen¬ 
dent  of  shift,  rotation,  size  and  representation.  They  are  the 
invariants  of  the  triangular  shape,  and  we  have  just  used  them 
to  define  equal  shape. 

A  check  of  the  literature  on  human  perception  shows  that 
experiments  exist  that  support  the  claimed  importance  of  the 
angles.  Ihe  eye  of  the  human  newborn  is  attracted  primarily 
by  the  angles  of  triangles  presented  to  it,  and  not  by  the 
sides  or  the  interior  (Salapatek  and  Kessen  1966).  A  further 
clue  that  the  eye  recognizes  patterns  from  their  contours  - 
and  particularly  from  the  change  of  the  contour  -  rather  than 
by  cross-correlation  with  sample  patterns  is  provided  by  the 
well-known  sensitivity  to  contours  and  the  insensitivity  to 
average  brightness.  A  contour  in  a  photograph  made  fuzzy  by 
motion  or  improper  focusing  is  spotted  immediately,  even  though 
it  would  have  almost  no  effect  on  the  integral  of  a  cross¬ 
correlation  with  a  sample  pattern.  On  the  other  hand,  the 
average  brightness  of  a  photograph  has  a  huge  effect  on  this 
integral,  but  the  eye  is  almost  insensitive  to  average  bright¬ 
ness  within  a  wide  range. 
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RECOGNITION  OF  POLYGONS 


Let  us  generalize  the  concept  of  invariants  of  shape. 

Figure  2a  shows  a  polygon  with  six  sides  a  ...  f  and  six  angles 
a  .  ..  £  .  Knowledge  of  the  angles  is  no  longer  sufficient,  we 
must  also  know  the  ratios  of  the  length  of  the  sides.  In 
order  to  arrive  at  a  unique  set  for  these  ratios,  we  look  first 
for  the  longest  side,  which  is  f.  We  then  form  the  ratios  a/f, 
b/f,  ...,  g/f,  1.  All  ratios  are  1  or  smaller.  Instead  of  the 
longest  side  one  may  also  use  the  sum  of  all  sides  for  normali¬ 
zation. 

The  longest  side  may  be  used  to  produce  a  unique  listing 
of  the  ratios  of  the  sides  as  well  as  of  the  angles.  We  start 
at  the  longest  side  f,  and  run  in  the  mathematically  positive 
sense  (counter  clockwise)  around  the  polygon.  The  following 
listing  of  side  ratios  and  angles  is  obtained: 

1,  a  ,  a/f,  P  ,  b/f,  . . e/f,  £  (1) 

Note  that  the  angles  a,  3,  ...  are  always  measured  in  the 
mathematically  positive  sense,  starting  from  the  side  already 
listed  and  ending  at  the  side  to  be  listed  next. 

If  two  polygons  produce  the  same  list,  we  say  they  have 
the  same  shape. 

A  couple  of  refinements  are  still,  needed.  First,  consider 
a  polygon  that  is  so  large  that  we  cannot  see  all  of  it  at  the 
same  time.  We  will  not  know  in  which  direction  to  proceed  in 
order  to  run  counter  clockwise  around  it.  A  machine  is  in  this 
position.  It  can  observe  only  local  features  of  a  polygon, 
either  a  straight  line  or  an  angle,  but  it  does  not  observe 
the  whole  polygon  at  once.  Without  knowing  what  will  turn  out 
to  be  running  around  the  polygon  in  the  mathematically  positive 
or  negative  sense,  one  has  no  choice  but  to  run  around  in  both 
possible  directions.  This  is  shown  in  Fig.  2b.  One  obtains 
the  listing 

1.  £  »  e/f,  e  ,  d/f,  6 . a/f, a  (2) 
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The  bars  over  the  angles  indicate  the  complement,  e.g.,  a  «Jnc. 

If  two  polygons  produce  the  same  listing,  thev  have  the 
same  shape.  Actually,  only  one  listing  has  to  be  produced  hv 
a  shape  recognizing  machine,  if  both  lists  are  stored  for 
re  ference . 


Let  us  replace  the  angles  in  Hqs.  (1)  and  (2)  by  their 
complement  s : 


l  ,  c.  ,  a/f,P,  b/  f ,  Y  »  •••,  e/f,£  (3) 

1 ,  £  ,  e/f  ,  e  ,  J/f  ,  6  ,  .  .  .  ,  a/f ,  a 


A  comparison  with  Figs.  2c  and  d  shows  that  these  two 
listings  define  the  mirror  image  of  the  shape  of  Fig.  2a.  Hence, 
complementing  the  angles  is  the  way  to  include  mirror  images 
into  the  definition  of  "equal  shape". 

The  .econd  required  refinement  is  for  polygons  that  have 
two  or  more  sides  of  equal  length,  and  which  are  at  the  same 
time  longer  than  all  the  other  sides.  An  example  is  shown  in 
Fig.  3,  where  the  two  sides  a  and  f  are  equally  long.  Starting 
with  f  and  proceeding  in  the  direction  indicated  in  Fig.  3a  we 
obtain  the  following  listing: 


1 ,  a  ,  1 ,  p  ,  b/ f ,  v ,  c/ f ,  6  ,  d  /  f ,  e 


e/f,  C 


14) 


However,  one  may  start  just  as  well  with  a  to  obtain  a  different 
list ing : 


l ,  P  ,  b/f,  y  ,  e/f,  6 


d/ f ,  *  ,  e/f,  ",  ,  1 ,  t 


IS) 


Since  F.qs.  (4)  and  (5)  differ  by  a  cyclic  shift  only,  one  has  to 
store  only  one  listing  for  reference.  If  the  machine  makes  a 
listing  of  a  polygon  that  has  to  be  identified,  and  more  than  one 
ratio  1  occurs  for  the  sides,  one  must  shift  the  listing  cycli¬ 
cally  so  that  each  l  comes  to  the  head  of  the  list,  and  compare 
each  such  list  with  the  stored  reference  list. 

Fig.  3b  repeats  Fig.  3a  for  running  around  the  polygon 
clockwise.  No  additional  problems  are  encountered. 
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When  the  shape  of  a  polygon  is  to  be  recognized  by  our  line 
extraction  method,  one  must  remember  the  starting  point,  since 
one  would  follow  otherwise  the  closed  line  forever.  Apart  from 
this  requirement,  our  method  may  be  used  in  exactly  the  same 
way  to  recognize  open  ended,  straight-line  shapes,  as  shown  in 
Figs.  4a  and  b.  However,  the  listing  cannot  start  with  the 


longest  line  but  with  one  of  the  lines  at  the  ends: 

a/b,  a  ,  1 ,  (3  ,  c/b  (6) 

c/b,  p  ,  1 , a  ,  a/b  (7) 

If  we  take  the  complements  of  the  angles,  we  obtain: 

a/b,  d ,  1,(3,  c/b  (8) 

c/b,  6  ,  l,  a  ,  a/b  (9) 


These  listings  define  the  mirror  image  shapes  of  Figs.  4c  and 
d. 

In  order  to  extend  the  theory  beyond  straight-line  struc¬ 
tures  we  need  a  generalization  of  the  concept  of  angle  that  is 
applicable  to  curved  lines.  Such  a  concept  is  the  curvature  of 
a  line.  The  curvature  is  invariant  to  shifts  and  rotations,  and 
by  introducing  a  normalization  factor  we  can  make  it  invariant 
to  size  or  scaling. 

Refer  to  Fig.  Sa  for  the  definition  of  curvature.  Two 
straight  lines  are  drawn  through  the  three  points  n-1,  n  and 
n*l.  The  angle  between  the  lines  at  the  point  n  is  denoted 
Ae(n),  and  the  distance  between  the  points  n  and  n+1  by  As(n). 
The  limit 

lim  At(n)/As(n)  =  dc/ds  (10) 

As(n)  — 0 


is  called  the  curvature  in  differential  mathematics. 


It  i  .  obvious  that  the  curvature  is  independent  of  any  shift 
or  rotation  of  the  three  points  n-1,  n  and  n*l.  Scaling  influence 
the  curvature  through  As(n),  but  a  multiplication  with  a  length 
characteristic  for  the  contour  will  remove  this  influence. 

The  main  difficulty  we  have  to  face  is  the  fact  that  an 
aerial  photograph  or  some  automated  analyzer  do  not  have  infinite 
resolution  and  thus  cannot  use  the  definition  of  Fq .  (10).  One 
mus i  develop  mathematical  methods  based  on  the  mathematics  of 
finite  differences.  To  show  why  this  must  lead  to  results  very 
different  from  the  ones  of  differential  mathematics,  we  inves¬ 
tigate  the  concept  of  "radius  of  curvature". 

In  fig.  5b  a  circle  with  radius  p(n)  is  drawn  through  the 
points  n  and  n+1.  The  line  from  n-1  to  n  equals  the  tangent 
of  the  circle  in  point  n.  We  can  determine  p(n)  and  define 
a  curvature  X(n)  =  1  /  P  ( n )  as  shown  in  fig.  5b. 

Instead  of  defining  a  circle  in  this  way,  we  may  also 
draw  a  circle  through  all  three  points  n-1,  n  and  n+1  as  shown 
in  Fig.  5c.  The  radius  p(n)  cannot  be  calculated  exactly  with¬ 
out  solving  a  transcendental  equation.  This  difficulty  does  not 
occur  in  differential  mathematics,  since  the  approximation 
sin  Ac  4 Ac  becomes  exact  when  Ac  becomes  a  differential  dc  . 
Furthermore,  in  differential  mathematics  we  would  have 

lim  s[As(n-l)  ♦  As(nl]  =  ds,  (11) 

As— >0 


and  the  relation 

p=  ds/dc  (12) 

would  follow  from  Fig.  5c,  which  is  just  the  inverse  of  Eq.(10) 
No  such  simple  relation  exists  for  finite  differences  As(n)  and 
Ac(n).  If  one  wants  to  define  a  radius  of  curvature  in  finite 
mathematics,  one  would  not  want  to  use  the  one  of  Fig.  5c  which 
requires  the  solution  of  a  transcendental  equation  but  prefer 
the  one  of  Fig.  5b  or  -  as  will  be  shown  later  on  -  a  still 


other  definition  that  is  more  useful  than  either  the  definition 
of  Fig.  5b  or  c. 


3.  CONTOUR  RECOGNITION 

We  base  the  following  investigation  on  the  premise  that 
a  pattern  or  a  contour  does  not  consist  of  nondenumcrably 
infinite  points  with  the  infinitesimal  area  dxdy  but  of  a  finite 
number  of  points  with  a  finite  area.  We  choose  this  finite 
area  to  be  squares  with  the  sides  Ax  and  Ay=  Ax,  having  thus  the 
area  AxAy-  An  ellipse  is  shown  in  this  representation  in  Fig.  6. 

The  reason  for  choos ing  squares  for  the  elemental  points  is 
that  liquid  crystal  displays  produce  such  points  (Harmuth  1977, 
Figs.  234-2  and  234-3).  Plasma  tube  displays  also  come  close 
to  producing  such  points.  A  cathode  ray  display  produces  cir¬ 
cular  rather  than  square  points,  and  these  points  are  not  in¬ 
herently  at  such  fixed  positions  as  in  liquid  crystal  and  plasma 
tvibe  displays. 

Let  the  34  x  26  elemental  points  in  Fig.  6  represent  a 
section  of  a  liauid  crystal  or  plasma  tube  display.  A  "scanning 
pattern”  of  four  elemental  points  -  as  shown  in  the  lower  left 
of  Fig.  6  -  shall  be  illuminated.  This  scanning  pattern  can 
be  moved  to  any  location  of  the  field  of  34  x  26  elemental 
points  in  Fig.  b.  In  the  case  of  a  liquid  crystal,  the  four 
points  of  the  scanning  pattern  arc  light  transparent  while  all 
the  other  points  of  the  field  of  34  x  26  points  are  opaque;  in 
the  case  of  a  plasma  tube  display  the  scanning  pattern  consists 
of  four  glow  discharge  points  while  all  the  other  points  are  dark 

Let  a  photo  transparency  be  overlayed  on  the  display. 

This  transparency  shall  be  opaque  for  the  approximately  ellip¬ 
tical  area  in  Fig.  6  containing  the  crossed  squares,  and  trans¬ 
parent  for  the  other  squares.  In  reality,  it  is  only  necessary 
that  the  contour  of  the  elliptical  area  is  noticeably  less 
transparent  than  the  area  outside  the  ellipse. 

The  next  step  is  to  make  the  scanning  pattern  run  around 
the  contour  of  the  ellipse  in  a  way  that  can  be  automated.  The 
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scanning  pattern  in  Fig.  6  is  moved  to  the  right  so  that  its 
point  0  coincides  with  the  point  1  =  0,  T)=l  of  the  display;  the 
scanning  pattern  is  then  in  the  location  £=0,  "0=1.  None  of 
the  four  illuminated  points  of  the  scanning  pattern  will  be 
blocked  by  an  opaque  section  of  the  transparency;  an  observer 
will  see  four  bright  points  or  the  pattern  denoted  0  in  the 
upper  left  corner  of  Fig.  7.  The  number  0  indicates  that  none 
of  the  2x2  illuminated  points  of  the  scanning  pattern  is 
blocked  by  an  opaque  section  of  the  transparency. 

The  pattern  0  in  Fig.  7  shows  the  command  "scan".  It 
means  that  the  scanning  pattern  should  be  moved  one  step  to 
the  right;  if  the  right  edge  of  the  field  in  Fig.  6  is  reached, 
"scan"  means  that  the  scanning  pattern  is  returned  to  the  left 
edge,  £  =  0,  and  moved  one  step  up  to  t|=2. 

Moving  the  scanning  pattern  one  step  to  the  right  from 
£  =  0,  7i  =  l  brings  it  to  the  point  £  =  1,  71  =  1.  Again,  the  pattern 
is  not  blocked  by  any  opaque  section,  the  contour  operator  0 
of  Fig.  7  and  the  command  "scan"  is  obtained.  This  is  repeated 
until  the  scanning  pattern  reaches  the  position  £  =  11,  T)=l; 
the  contour  operator  IB  of  Fig.  7  is  now  obtained. 

Contour  operator  IB  is  shown  with  Arrow  1  pointing  to  the 
right  and  Arrow  2  pointing  upward.  This  means  that  the  scanning 
pattern  should  as  the  first  choice  be  moved  +1  in  the  direction 
of  £.  If  this  first  choice  is  prohibited,  the  scanning  pattern 
is  moved  as  the  second  choice  +1  in  the  direction  of  Tl.  The 
only  reason  for  not  using  the  first  choice  is  that  the  scanning 
pattern  would  be  moved  back  to  the  point  it  had  just  previously 
been.  A  practical  example  will  show  the  working  of  the  contour 
operators  better  than  many  words. 

Table  1  lists  on  the  top  left  the  valiies  £=11,  7]=1  for 
which  the  contour  operator  IB  was  observed.  The  point  n=l  is 
recognized,  and  Fig.  7  shows  that  one  has  to  advance  +1  in  the 
direction  £  if  the  contour  operator  IB  is  observed.  This  brings  the 
scanning  pattern  to  £=12,  Tl  =  l.  The  contour  operator  2B  of 
Fig.  7  is  observed,  the  point  n=2  is  recognized,  and  one  must 
advance  +1  in  direction  £.  The  advancement  by  +1  in  direction 
£  and  the  contour  operator  2B  are  repeated  until  the  point 
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Table  1.  Scanning  of  the  contour  of  the  ellipse  ol' 


£  *21 ,  ')  •  1  is  reached .  Now  the  contour  operator  111  is  observed, 
no  new  point  is  recognized,  and  the  scanning  pattern  must  be 
advanced  M  in  the  direction  * .  We  leave  it  to  the  reader  to 
verify  the  entries  in  the  first  six  columns  of  Table  1  with  the 
help  of  Figs,  b  and  7. 

The  entries  (. ♦  1  )  in  the  column  "command  to  advance,  ti  "  in 
Table  1  signify  that  the  second  choice  rather  than  the  first 
choice  of  advancement  is  used,  which  would  be  to  advance  *1 
in  direction  F,  according  to  Fig.  7.  One  may  readily  verify  that 
this  first  choice  would  return  the  scanning  pattern  to  the 
point  it  had  just  been  previously  (from  h,x29,  ti=l5  to  {j-30, 

“IS;  from  1,  "28,  '1  - 1  7  to  1,-29,  11-17;  etc.). 

The  three  additional  columns  of  Table  I  will  be  discussed 
later  on  in  more  detail.  They  show  the  angle  Ay(n)  between 
two  lines  drawn  through  the  points  n  1,  n  and  n,n+l,  furthermore 
the  line  element  As  (n  )/Ax  and  a  curvature  Ay(n)Ax/A s  (n) . 

Let  us  discuss  how  a  machine  can  recognize  the  contour 
operators  of  Fig.  7.  In  principle,  one  could  do  just  like  our 
eyes  do  and  use  many  -  or  at  least  four  -  receptors  to  recognize 
the  operation.  However,  this  leads  to  great  problems  in  the 
implementation  of  equipment.  A  better  solution  is  to  make  a 
sequence  of  four  detections  rather  than  four  simultaneous  detec¬ 
tions.  The  reason  is  not  so  much  to  reduce  the  number  of  photo 
detectors  and  amplifying  circuits  from  four  to  one,  but  to 
simplify  the  optical  design. 

Figure  8  shows  again  the  scanning  pattern  with  four  bright 
points  of  Fig.  6  in  the  upper  left  corner.  This  pattern  is 
now  denoted  as  the  two-dimensional  Walsh  function  wal(0,|)X 
wa  1  (  0 ,  n  ) .  In  addition,  the  Walsh  functions  wal(0,|)  wa  1  ( 1  ,T|  ) , 
wa  l  (1,1.)  wa  1  (. 0 ,  rj )  and  wa  1(1,1, )  wa  1  (1,ij)  are  shown.  The  white 
areas  represent  the  value  M,  the  black  areas  the  value  -1.  We 
produce  the  Walsh  transform  of  the  lb  contour  operators  with 
these  four  functions.  In  terms  of  differential  mathematics  we 
get  the  coefficients  a(i,k) 

■  I  7  F(l,  ,i|  )  wa  I  (i  ,1,  )  wal  ( k .  r) )  dl,dii 
0  0 


-  in  - 


a  (  i  .  k  ) 


(13) 


We  have  to  insert  each  one  of  the  16  contour  operators  of  Fig. 

7  foT  F(|,ti).  Ideally,  the  crossed  squares  of  the  operator  are 
perfectly  opaque  and  their  brightness  equals  0;  the  clear  squares 
are  fully  transparent  and  their  brightness  may  be  set  equal  to 
+1.  In  reality,  we  would  like  to  recognize  the  elliptical  area 
in  Fig.  6  if  it  is  somewhat  less  transparent  than  the  surround¬ 
ing  area,  since  this  is  what  our  eyes  can  accomplish.  Hence,  we 
assign  the  brightness  q  to  the  crossed,  opaque  squares  of  the 
operators  in  Fig.  7,  and  the  brightness  p  to  the  clear  squares: 

p  =  brightness  of  transparent  elemental  point  (14) 


q  =  brightness  of  opaque  or  less  transparent  elemental  point 
p/q  =  r  -  brightness  ratio 


0  =  q  <  1 ,  0  <  p  =  1 ,  p/q  >  1 


In  terms  of  finite  mathematics  we  represent  each  operator  of 
Fig.  7  by  an  array  of  four  coefficients  p  and  q;  e.g.,  the 
operator  1A  is  represented  as  follows 


P  P 
P  q 


(15) 


The  four  Walsh  functions  of  Fig.  8  are  represented  by  four  arrays 
with  the  coefficients  +1  and  -1: 


■+1  +1 

■-i  -r 

+1  +1 

+i  +i 

+i  -i 

"-i  +i 

+i  -i 

+i  -i 

(16) 


Multiplication  of  the  array  of  Eq.(15)  with  the  arrays  of  Eq.(16) 
and  summation  of  the  products  yield  the  coefficient  a(i,k): 


a(0,0) 


+  1 
+  1 


3p  +  q  =  3p (1+q/p) 


(17) 
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a(0,l)  -V 

P  P 

p  q 

X 

-1 

+1 

-l1 

+1 

=  -p+  q  =  -p  (1  -q/p) 

ad.o)  =>; 

P  P 

X 

+  1 

-1 

=  p  -  q  =  p (1 -q/p) 

a  (1,1) 

P  <1. 
"p  P 

X 

+1 

■-1 

-1 

=  p  -  q  =  p (1 -q/p) 

_p  q. 

_+l  -1J 

The  multiplication  and  summation  process  is  the  same  as  used  for 
the  reduction  of  two  tensors  of  rank  2  to  a  constant. 

The  four  coefficients  of  Eq.(17)  are  listed  in  Table  2  for 
the  operator  1A.  The  listed  coefficients  for  the  other  operators 
are  obtained  in  the  same  way. 

Since  the  absolute  brightness  of  p  and  q  is  of  less  interest 
than  their  ratio  r  =  p/q,  we  divide  the  coefficients  a(0,l), 
a(l,0)  and  a(l,l)  by  a(0,0).  The  values  in  columns  6-8  of 
Table  2  are  obtained.  One  may  verify  that  R/3  is  always  smaller 
than  Q  or  P.  Hence,  by  setting  threshholds  at  +R/6  and  -R/6 
one  can  divide  the  coefficients  in  columns  6-8  into  three  classes: 


a)  coefficients  larger  than  +R/6 

b)  coefficients  between  +R/6  and  -R/6 

c)  coefficients  smaller  than  -R/6 


(18) 


The  numbers  +  1,0  and  '1  are  written  in  columns  9-11  of  Table  2 
for  these  three  classes  of  coefficients.  Each  set  of  these 
"quantized"  coefficients  defines  a  contour  operator,  e.g.,  the 
set  -1,0,0  defines  the  operator  2D.  Hence,  by  using  the  four 
Walsh  patterns  of  Fig.  8a  and  doing  the  processing  according  to 
Table  2  we  need  only  one  phototube  with  associated  circuits  but 
four  observations  to  recognize  any  one  of  the  16  contour  operators 
in  Fig.  7  uniquely. 

Since  a  liquid  crystal  or  a  plasma  tube  produces  either 
"light"  or  "no  light"  at  a  point,  but  never  "positive  light"  and 
"negative  light",  we  still  have  to  explain  how  we  can  implement 
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the  Walsh  functions  practically.  Actually,  we  do  not  need  these 
functions  but  the  coefficients  a(0,0),  a(0,l),  a(l,0),  and  a(l,l) 
of  Table  2.  We  may  use  the  two-dimensional  block  functions 

bio  ( i ,  ^ )  bio  (k,t))  of  Fig.  8b  to  obtain  the  Walsh  coefficients. 

v-/ 

Let  the  elemental  point  represented  by  the  clear  square  of 
the  function  bio  (0,!-)  bio  ( 0 ,  T| )  be  illuminated  and  the  three 
elemental  points  represented  by  the  black  squares  be  dark.  This 
implies  that  one  square  of  a  liquid  crystal  is  transparent  or 
one  glow  point  of  a  plasma  tube  emits  light.  A  photo  transparency 
be  overlayed;  one  of  the  sixteen  patterns  of  Fig.  7  is  thus  over- 
layed  on  the  pattern  bio  (0,£)  bio  ( 0 ,  T) )  of  Fig.  8.  The  light 
passing  through  to  the  phototube  will  be  proportionate  to  either 
p  or  q,  depending  on  whether  a  transparent  or  an  opaque  square 
of  the  photo  transparency  overlays  the  bright  square  of  bio  (0,^)x 
bio  (0,t|).  Whichever  it  is,  the  output  voltage  of  the  phototube 
represents  the  coefficient  b(0,0),  which  can  be  written  in 
analogy  to  Eq.  (13)  : 


7  2r  -  ^ 

b  ( i ,  k)  «  J  £  F(£,ti)  bio  (i,|)  bio  (k,Ti)  d^dt) 


(19) 


Illumination  of  the  liquid  crystal  or  the  plasma  tube  according 
to  the  other  three  functions  of  Fig.  8b  produces  the  coefficients 
b(0,l),  b(l,0)  and  b ( 1 , 1 ) .  A  comparison  of  the  Walsh  and  the 
block  functions  yields  immediately  the  following  relations: 


a (0 , 0)  »  b (0 , 0)  ♦  b ( 0 , 1 )  ♦  b (1 , 0)  +  b ( 1 , 1 ) 


(20) 


a (0 , 1 )  =  b (0 , 0)  ♦  b(0,l)  -  b(l,0)  -  b(l,l) 


a  (1 , 0)  -  b (0 , 0)  -  b (0 , 1 )  ♦  b(l,0)  -  b(l , 1) 


a(l  ,1)  *  b ( 0 , 0 )  -  b(0, 1)  -  b (1 , 0)  ♦  b(l,l) 


-  14  - 


Hence,  we  have  a  practical  way  to  obtain  the  coefficients 
a(0,0)  to  a(l,l)  of  Table  2.  There  is  another  way  that  does 
not  require  the  block  functions;  it  is  based  on  interpreting 
the  black  areas  ol  the  Walsh  functions  as  "no  light  emitted"  and 
the  clear  areas  as  "light  emitted." 

Summing  up,  we  have  developed  a  deterministic  method  to 
find  point  -  by-point  the  contour  between  areas  of  different 
light  transmission  of  a  photo  transparency,  if  this  transparency 
is  placed  on  a  liquid  crystal  or  a  plasma  tube  display.  The 
method  requires  that  one  point  of  the  display  emits  light, 
which  passes  through  the  transparency  to  the  phototube.  A 
sequence  of  four  signals  received  by  the  phototube  characterizes 
one  of  the  16  contour  operators  of  Fig.  7,  and  this  operator 
commands  which  quadruple  of  points  should  be  illuminated  next. 

Let  us  end  this  section  with  an  explanation  how  the  arrows 
for  the  advancement  were  chosen  in  Fig.  7.  The  two  arrows  shown 
for  the  contour  operator  1A  assure  that  the  opaque  square  will 
occur  in  the  following  operator  too.  For  instance,  if  the 
operator  1A  is  recognized  and  the  scanning  pattern  is  advanced 
one  step  to  the  right,  the  next  recognized  pattern  could  be 
either  IP,  2D,  2E  or  3D.  If,  on  the  other  hand,  the  scanning 
pattern  would  have  advanced  one  step  to  the  left  after  recog¬ 
nition  of  the  contour  operator  IA,  one  might  have  lost  the 
opaque  point  and  obtained  the  contour  pattern  0,  which  means 
the  contour  would  have  been  lost  and  we  would  be  back  to  scan¬ 
ning  to  find  a  contour. 

One  may  readily  verify  that  all  arrows  in  Fig.  7  have  been 
chosen  so  that  the  scanning  pattern  will  never  advance  into  an 
area  that  yields  either  the  all  transparent  operator  0  or  the  all 
opaque  operator  4.  Putting  it  differently,  the  arrows  arc 
chosen  so  that  the  next  contour  operator  will  have  at  least  one 
transparent  square  and  at  least  one  opaque  square. 

The  choice  of  the  advance  in  the  directions  of  the  two 
arrows  is  arbitrary.  We  have  chosen  the  following  order  of 
preference:  1.  advance  right;  2.  advance  up;  3.  advance 

left;  4.  advance  down.  In  other  words,  the  arrow  having  the 
smaller  angle  with  the  positive  x-axis  got  the  number  1,  the 


other  got  the  number  2.  This  rule  also  applies  to  the  operators 
2E  and  2F. 

4.  CURVATURE  IN  FINITE  MATHEMATICS 


We  have  seen  in  Fig.  5a  that  one  needs  three  adjacent  points 
to  define  the  angle  AE(n)  and  the  curvature  in  point  n.  The 
three  points  n-1,  n  and  n+1  could  have  any  position  relative 
to  each  other  in  Fig.  5,  but  the  square  pattern  of  Fig.  6  res¬ 
trains  this  freedom  greatly.  Figure  9  shows  three  typical 
examples  lor  the  location  of  three  adjacent  points  in  the  square 
field  of  I  ig.  o.  We  must  develop  some  concepts  of  finite  mathe¬ 
matics  on  a  square  field  before  we  can  proceed. 

In  differential  mathematics  using  cartesian  coordinates  the 
curvature  is  defined  by  the  relation 


xlx) 


dr 

Ts 


(W2)577 


(a) 


This  expression  requires  the  first  and  second  derivatives,  and 
we  will  thus  investigate  the  first  and  second  derivatives  of 
finite  mathematics  on  a  square  field. 

In  finite  mathematics  we  have  at  the  point  n  of  a  curve 
a  left  difference  quotient  t^(n)  and  a  right  difference  quotient 
Ti;(n). 


^(n) 


1)1 
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i  -  ill 

Ln- 1) 

51 

ini 

1  -  H 

[TT-TT 

Hj.00 


(n*l] 

1  -  ill 

S''}} 

(n+1  ] 

'  - 

n- 1 ) 

122) 

12  5) 


Refer  to  Fig.  9a.  The  variable  n  determines  a  point  along  the 
line  by  counting  from  the  arbitrarily  chosen  first  point  n«l  . 
The  equivalent  in  differential  mathematics  is  the  length  of  the 
curve  from  the  arbitrarily  chosen  beginning  to  the  point  n. 

The  cartesian  coordinates  of  this  point  n  are  £ -j  and 
The  left  difference  quotient  becomes 


-  16  - 


-  0, 


(k-1)  -  lk-1) 

j  -  n -n 


and  the  right  difference  quotient  becomes 


^(n)  -  nvrp---1}-  m  +1 

The  two  quotients  and  T)j.  are  shown  in  Fig.  9a  to  determine 
the  tangent  of  the  angles  y^(.n)  and  vr(n)  with  which  straight 
lines  through  the  points  n-1,  n  and  n,  n+1  intersect  the 
j'-axis.  The  angle  between  the  two  lines  is  denoted  A>(n): 

A  V  (n )  -  y  r  (n)  -  y£(n)  (-1) 

Figures  9b  and  c  give  two  more  examples  of  three  adjacent 
points,  and  the  two  traight  lines  with  the  angle  Ay  between 
them  defined  by  three  points.  A  closer  study  shows  that  there 
are  only  40  possible  sequences  of  three  points  on  a  square 
field.  The  first  twenty  are  shown  in  Fig.  10,  the  other  twenty 
in  Fig.  11.  These  two  figures  differ  only  in  the  direction  in 
which  the  numbering  of  the  three  adjacent  points  increases.  A 
;ero  ”0"  i  r.  shown  to  indicate  this  direction  of  numbering. 

Only  eight  different  values  of  the  left  and  the  right 
difference  quotients  ^  and  T|^  occur  in  Figs.  10  and  11.  These 
values  are  listed  in  Tables  A  and  4.  A  complete  listing  of  the 
quotients  and  t] '  for  the  40  configurations  a  -  t  and  a'  *  t' 
of  Figs.  10  and  11  is  given  in  Table  5.  Also  shown  is  the  second 
order  difference  quotient  t]"ln): 


fl "  (.  n ) 


’1 

t 


C4) 


Finally  the  angle  Ay(n)  defined  by  the  two  straight  lines  through 
the  three  points  -  intersecting  in  the  center  point  -  are  shown. 
One  may  see  that  r\£ ,  T)^  and  t)H  uniquely  define  Ay.  Since  the 
three  difference  quotients  in  turn  are  uniquely  defined  by 
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T  t  le  l- .  The  eight  possible  values  of  the  right  difference  quotient. 


'rv  J^n+1  ^-guO 


I(n-l),  £(n),  £(n+l),  m(n-l),  i](n)  and  r)(n*l),  which  have  been 
obtained  in  Table  1  by  means  of  the  contour  operators,  we  can 
derive  the  values  of  by  means  of  the  contour  operators. 

Let  us  show  an  example.  For  point  n-10  in  column  4  of 
Table  1  we  find  the  values  £(9)  =  19,  £(.10)  =  20,  £(11)  =  21, 

Tl(9)  «  1,  >1(10)  =  l  and  T](  11)  =  2.  It  follows: 

"i1101  ’  .4^  '  »•  ’ijuo)  ■  ■  i.  i"oo)  -  frrrjvrs 

Table  5  shows  this  set  of  values  for  the  sequence  e  in  Fig.  10, 
and  the  angle  A  y  =  Tt/4  in  correspondence  with  Table  1. 

As  a  second  example,  consider  the  point  n=ll  of  Table  1  with 
£(10)  =  20,  CUD  =  21,  £(12)  =  22,  q(l  0)  =  1,  ri(ll)  =  2  and 
11(12)  =■  2. 

i^UD  -  yfry^  =  U  ii’UD  -  zl-Tl  =  °»  ^"UU  =  yrril^Io  ■50° 

This  set  of  values  holds  for  the  sequnce  h  in  Table  5  and  Fig.  10. 
The  angle  AY  =  -it/ 4  is  obtained  in  agreement  with  Ay(ll)  in 
Table  1. 

The  definition  of  Eq.(10)  for  the  curvature  contained  the 

angle  A  £ ( n )  -  which  we  denote  now  by  Ay(n)  for  reasons  to  be 

seen  later  -  and  the  length  element  As(n).  We  will  not  use  the 

length  clement  in  this  form  but  for  completeness  sake  we  show  in 

Fig.  12  that  in  finite  mathematics  on  a  square  field  there  can 

be  only  three  values  for  the  normalized  length  element  As(n)/Ax. 

The  value  of  As(n)/Ax  is  uniquely  defined  by  the  quotients  rj^(n), 

Ti'(n)  and  V'(n),  since  one  could  list  the  value  of  As(n)/Ax  in 
r  1 

addition  to  the  value  of  Ay(n)  in  Table  5. 

The  straight  forward  way  to  arrive  at  a  curvature  as  function 
of  the  point  n  would  be  to  calculate  the  line  element  As(n)/Ax  and 
the  curvature  Ay (n)Ax/A s (n)  as  shown  in  the  last  two  columns  of 
Table  1.  It  can  be  seen  that  this  is  not  a  satisfactory  way 
since  the  curvature  in  Table  1  is  mostly  0  and  alternates  between 
positive  and  negative  values.  One  would  have  to  average  the 
curvature  over  many  values  of  n  to  obtain  a  satisfactory  result. 
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The  first  reason  for  not  following  the  straight  forward  way  is 
thus  that  averaging  is  an  unsatisfactory  method.  There  is, 
however,  an  even  more  important  reason.  In  aerial  photography 
analysis  we  are  mainly  interested  in  man-made  objects.  The 
contours  of  such  objects  are  primarily  straight  lines  and  circles. 
Hence,  one  will  want  a  method  that  first  tries  to  fit  the  contour 
obtained  from  the  contour  operators  by  straight  lines;  if  this 
is  impossible,  a  fit  by  circles  should  be  tried  next;  and  only 
if  straight  lines  and  circles  fail  will  one  admit  a  more  general 
curvature . 

For  an  illustration  of  the  importance  of  this  concept  refer 
to  Fig.  13,  which  shows  an  opaque  rectangular  area  on  the  same 
square  field  as  the  ellipse  of  Fig.  6.  Table  6  shows  the  charac¬ 
teristic  values  for  the  scanning  of  the  contour  of  the  rectangle 
in  analogy  to  Table  1.  I f  we  simply  average  the  values  for  the 
curvature  in  the  last  column,  we  would  obtain  zero  curvature  for 
the  straight  lines  of  the  rectangle,  but  the  corners  at  points 
n=l,  26,  39,  64  would  be  rounded. 

5.  CONTINUOUS  CURVE,  DISCONTINUOUS  FIRST  DERIVATIVE 

We  return  to  Table  1  but  we  disregard  the  values  listed  for 
the  line  element  and  the  curvature.  Using  only  the  values  of 
A{-  and  Ati  as  function  of  the  variable  n  we  can  plot  the  72  points 
n=l  to  n=72  of  Fig.  14a.  This  plot  is  already  independent  of  any 
shifts  in  the  directions  and  T]  as  well  as  from  rotations.  The 
independence  from  shifts  is  evident  from  the  way  Al,  and  At]  were 
obtained,  while  the  independence  from  rotations  is  not  so  clear. 

We  will  demonstrate  this  independence  later  on,  and  it  will 
become  clear  that  the  independence  is  obscured  by  the  poor  resolu¬ 
tion  of  our  ellipse.  The  measure  for  resolution  is  the  number 
of  points  of  the  contour,  which  is  72  in  Fig.  14a.  For  comparison, 
the  eye  can  resolve  about  0.1  mm.  Hence,  75  points  would  imply 
an  ellipse  with  a  circumference  of  7.5  mm.  The  eye  would  have 
trouble  to  decide  whether  such  a  small  structure  is  an  ellipse. 

As  a  result  of  this  small  number  of  points,  our  results  will  not 
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Table  6.  Scanning  the  contour  or  the  rectangle  or  F ig.13. 
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look  very  good,  but  the  only  way  to  improve  them  is  to  increase 
the  resolution  by  increasing  the  number  of  points  by  an  order 
of  magnitude.  This  is  no  problem  if  the  methods  derived  here 
are  implemented  by  machines,  but  our  drawings  and  numerical 
examples  would  become  unmanageable  with  such  a  large  number  of 
points . 

At  this  time,  we  want  to  turn  to  the  problem  of  size  or 
scaling.  Consider  an  enlargement  or  a  reduction  of  the  "ellipse" 
represented  by  the  72  points  in  Fig.  14a.  What  would  be  enlarged 
or  reduced  would  be  the  distances  between  the  points,  but  their 
number  would  remain  unchanged.  On  the  other  hand,  it  is  quite 
evident  from  Fig.  6  that  a  larger  ellipse  would  yejild  more 
points  rather  than  larger  distances  between  the  same  number  of 
points.  Such  a  problem  does  not  occur  in  differential  mathe¬ 
matics  since  all  curves  have  there  nondenumerably  many  points, 
and  as  a  result  there  are  no  gaps  or  distances  between  the 
points . 

In  order  to  make  Fig.  14a  scalable,  we  have  to  replace  the 
finite  set  of  points  by  a  nondenumerable  set.  This  is  readily 
done  by  connecting  adjacent  points  by  straight  lines  as  shown 
in  Fig.  14.  Let  us  choose  the  coordinates  of  a  point  n  arbitra¬ 
rily  to  be  £  q,  tip.  The  coordinates  of  the  point  n+1  are  then 
+Ag ,  +  At)  .  The  values  A|  and  At)  come  for  any  pair  of 

points  n,  n+1  from  Tables  1  or  6.  The  equation  of  the  line 
between  the  two  points  is 

(E'Eo>  *  \  (25) 

Using  this  equation  with  arbitrary  initial  values  for  the 

point  n  and  arbitrary  direction  of  the  |-axis  -  the  T|-axis 
being  perpendicular  and  positive  oriented  -  a  computer  can  plot 
Fig.  14b  using  the  values  of  n,  A£  and  At)  in  Table  1  . 

Figure  14b  can  be  scaled  in  the  usual  meaning  of  this  term 
since  the  gaps  between  the  points  have  been  filled  by  lines  with 
nondenumerably  many  points.  We  have  now  a  continuous  curve 
with  a  discontinuous  first  derivative. 


A  human  observer  looking  at  Fig.  14b  and  knowing  that  it  rep¬ 
resents  a  contour  on  a  square  field  will  immediately  smooth  out 
the  indentations  at  the  points  n-1,  14,  20,  ...  to  obtain  the  smoothe 
contour  of  Fig.  15.  This  is  still  a  continuous  curve  with  discon¬ 
tinuous  first  derivative,  but  it  is  definitely  better  than  Fig.  14b. 
Mathematically  speaking,  Fig.  15  is  obtained  by  turning  a  straight 
line  around  every  point  in  a  positive ‘sense.  How  can  one  automate 
this  process  with  the  values  in  Table  1? 

We  first  rewrite  Table  1  in  a  concentrated  form.  We  do  no  longer 
need  the  "location  of  the  scanning  pattern"  and  the  "observed  contour 
point".  We  need  n,  A  F  ,  An  and  Ay.  Since  the  ellipse  of  Fig.  6 
contains  72  points,  the  values  of  AF,  At}  and  Ay  repeat  after  point 
36.  Hence,  we  list  A^fA"n  and  Ay/ 2rc ;  n=l  ...  36  in  Table  7.  Further 
more,  we  list  in  column  5  the  sum  of  Ay/2n  taken  from  n=l  to  n. 

We  see  that  this  sum  increases  at  the  points  n=nc  listed  in  column 
6.  A  comparison  with  Fig.  15  shows  that  these  are  just  the  points 
needed  for  the  construction  of  this  curve.  One  may,  of  course, 
just  as  well  say  that  the  points  nc  are  the  one  distinguished  by 
the  positive  values  ofAy/2rt  in  column  4. 

Next  we  need  the  angles  between  the  straight  lines  at  the 
points  n=nc  in  Fig.  15.  Refer  to  Fig.  16.  On  top  is  shown  a 
typical  contour  of  an  opaque  area  on  a  square  field.  The  polygon 
approximation  we  want  is  shown  by  the  dashed  line  through  the 
points  n=l,  2,  5,  7,  10,  11,  12.  Figure  16b  shows  the  contour 
points  alone.  The  dashed-dotted  line  connecting  all  points  corres¬ 
ponds  to  the  approximation  of  Fig.  14b.  The  angles  at  the  points 
n  are  denoted  Ay(n)  and  their  curve 

Y(n)  *  2  Ay (5.)  (26) 

i  =  l 

is  plotted  in  Fig.  16c.  The  angles  of  the  dashed  lines  with  a 

certain  reference  line  -  typically  the  F-axis  -  in  Fig.  16b  are 

denoted  e(n)  =  e(nc) .  The  difference  of  the  angles  in  a  point  n 

is  denoted  Ae(n).  If  we  denote  a  certain  corner  point  nc  with 

n  ,  and  the  next  larger  one  with  n  ,  we  obtain 
cl  c2 

Ar(nci)  -  e(nc2)  *  G^ncl)  (27) 

This  relationship  is  quite  evident  in  Fig.  16b  for  ncl  =  5  and 
nc2  =  7.  The  positive  jumps  of  y(n)  at  the  corner  points 
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nc  =  2,  5,  7,  10  is  also  evident. 

The  value  of  e(n)  follows  directly  from  an  inspection  of 
Fig.  17  as 

£(n)  ■  arc  tan  Y/X  (28) 

X  is  the  number  of  steps  A|  in  the  direction  |  and  Y  is  the 

number  of  steps  Atj  in  the  direction  7)  from  the  corner  point 

n=n  =n  .  to  the  following  corner  point  n=n  The  sums  X  and 
C  Cl  c  z 

Y  of  and  At)  are  listed  in  columns  6  and  7  of  Table  7  for  the 
ellipse  of  Fig.  6. 

We  see  further  from  Fig.  17  that  the  distance  As(nc)/Ax 
from  one  corner  point  to  the  next  equals 

As(nc)/Ax  =  (X2  +  y2)1/Z  (29) 

To  simplify  the  notation  we  will  generally  not  write  the  dis¬ 
tance  unit  Ax,  and  thus  simply  write  As(n)  for  As(n)/ 

The  values  of  As(n)/Ax  =As,e(n)  and  Ae(n)  are  listed  in 
columns  9-11  of  Table  7. 

We  still  need  the  total  length  of  the  contour  of  Fig.  15 
from  the  point  n=l  to  any  corner  point  n-nc  =  10,  13,  ....  It 
follows  from  Eq. (29) 

n*n  - 1 

s(nc)  =  2  As(i)  (30) 

i  =  1 

nc  =  1 ,  10,  13,  15  .... 

For  nc=l  the  upper  limit  of  this  sum  equals  zero  and  there 
is  thus  no  term  As(i)  in  the  sum;  hence,  s(l)=0.  For  nc=10 
we  have  the  upper  limit  n=9,  and  the  only  term  in  the  sum  is 
As(l)  =  9,  which  yields  s(10)  =  9.  The  values  of  s(nc)  are 
listed  in  column  12  of  Table  7.  The  largest  value,  s(37)  * 
40.58,  is  half  the  circumference  of  the  approximated  ellipse. 

In  order  to  scale  all  possible  ellipses  to  the  same 
reference  size,  we  multiply  their  circumference  with  a  factor 
that  makes  it  100  or  rather  lOOAx.  The  number  100  is  arbitrary, 
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and  it  is  chosen  because  it  is  small  enough  to  he  manageable. 

For  machine  operation,  one  would  choose  a  much  larger  number 
such  as  1000  or  10000.  It  should  be  as  large  or  larger  as  the 
largest  number  of  points  of  a  contour  in  a  square  field  such 
as  Figs.  6  and  13,  in  order  to  avoid  a  loss  of  information  by 
sea  1 ing . 

The  normalized  length  s*(n^.)  between  the  corner  points 
follows  from  the  relation 

s*(nc)  =  (50/40. S8)  s(nc)  (31) 

and  it  is  listed  in  column  13  of  Table  7. 

6.  CONTINUOUS  FIRST  DERIVATIVE,  DISCONTINUOUS  SECOND 

DERIVATIVE 

We  have  developed  enough  mathematics  for  the  automated 
plotting  of  the  approximation  of  an  ellipse  according  to  Fig.  15 
from  the  observed  contour  operators  of  Table  1.  We  could  now 
start  discussion  how  such  a  plot  could  be  recognized  by  a  com¬ 
puter  in  the  sense  that  the  printout  "ellipse  with  axis  ratio 
p  and  circumference  S  recognized"  is  produced  rather  than  the 
plot  of  Fig.  15.  However,  we  will  not  do  so  yet  for  the  follow¬ 
ing  reason.  We  know  from  experience  that  the  eye  is  very  sen¬ 
sitive  to  breaks  in  a  straight  line,  which  means  the  angles  in 
Fig.  15.  Mathematically  speaking,  a  discontinuity  of  the 
first  derivative  is  readily  spotted.  This  is  not  all.  The  eye 
is  also  sensitive  to  a  discontinuity  of  the  curvature.  This  is 
an  experience  gained  by  carefully  drafting  curves  for  publica¬ 
tion.  A  jump  of  the  tangent  in  a  point,  which  is  a  discontinuity 
of  the  first  derivative,  is  a  glaring  deficiency.  A  jump  of 
the  curvature  obtained,  e.g.,  by  joining  two  circular  sections 
with  different  radius  but  the  same  tangent  at  the  junction, is 
less  glaring  but  still  noticeable.  Hence,  we  will  carry  the 
theory  from  the  continuous  curve  with  discontinuous  first  deri¬ 
vative  to  a  continuous  curve  with  continuous  first  and  second 
derivative  before  we  turn  to  the  recognition  and  classification 
of  the  curve.  According  to  F.q.(21)  one  needs  a  continuous  first 
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and  second  derivative  to  obtain  a  continuous  curvature. 

Let  us  observe  at  this  time  that  the  straight  line  approx¬ 
imation  of  Fig.  15  is  already  sufficient  to  recognize  all 
contours  consisting  of  straight  lines,  such  as  the  rectangle 
in  Fig.  13.  To  show  this  more  clearly,  we  write  Table  8  for 
the  rectangle  of  Fig.  13  in  analogy  to  Table  7  for  the  ellipse 
of  Fig.  6.  The  interesting  values  are  Ae  and  s(nc).  With  them, 
derived  for  all  points  n=l  ...  76  rather  than  only  the  part 
shown  in  Table  8,  one  can  construct  Fig.  18.  The  rotation  of 
Fig.  13  is  removed,  and  there  are  clearly  four  straight  lines. 

The  corners  do  not  show  one  angle  Ae=Tt/2  but  either  two  or 
three  angles  that  sum  to  ti/2: 

Ae(26)  =  l .  2490,  At  (27)  =  0. 3218  ;  Ae(26)+Ae(27)*rt/2 

A  e(76 )  =  0 . 4636 ,  A e(  1 )  =  0 . 78 54  ,  Ae(2)«0. 3218  ;  Ae (76)  +  Ae(l)  +  Ae (2)  =ft/2 

There  is  no  particular  problem  for  a  machine  to  recognize  a 
rectangle  from  the  numbers  in  columns  11  and  12  of  Table  8, 
particularly  if  the  circumference  of  the  rectangle  contains  an 
order  of  magnitude  more  points  than  the  one  in  Fig.  13. 

Let  us  return  to  the  approximation  of  the  ellipse  in  Fig. 

15.  If  we  used  an  order  of  magnitude  more  points  we  would 
increase  the  number  of  corner  points  where  At  is  not  zero,  in 
contrast  to  the  rectangle  of  Fig.  18.  We  would  take  this  large 
number  of  corner  points  as  a  reason  to  decide  that  the  contour 
does  not  consist  of  straight  lines  but  is  curved.  Since  man¬ 
made  curved  contours  are  primarily  circles,  we  will  use  sec¬ 
tions  of  circles  to  smooth  the  corners  in  Fig.  15.  There  are 
many  ways  of  doing  so;  we  chose  the  one  shown  in  Fig.  19. 

The  three  points  n-1,  n  and  n  +  1  are  corner  points  in  our 
terminology.  The  distance  As(n)  is  smaller  than  the  distance 
As(n-l).  A  circle  is  drawn  from  the  center  of  the  line  between 

n  and  n+1  to  the  point  on  the  line  from  n-1  to  n  having  the 

distance  \s(n)/2  from  n.  The  lines  from  n  to  n+1  and  from 
n-1  to  n  are  the  tangents  of  the  circle  at  the  points  where 
the  lines  meet  the  circle.  This  insures  a  continuous  first 
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Table  8.  Computation  of  the  approximation  by  a  continuous  curve  with 
po.  itive  turning  tangent  and  discontinuous  first  derivative  for  the 


n 

A* 

73 

70 

►  1 

►  1  8 
►  1/8 
-  1/8 

0 

►  1/8 
-1/8 
0 

♦  1/8 
-1/8 

0 

►  1/8 
-1/8 

0 

►  1/8 
-1/8 
0 

♦  1  /8 
-1/8 
0 

♦  1/8 

-1/8 

0 

►  1/8 
-1/8 
0 

♦  1/4 

♦  1  /  8 
-  1/8 
0 

t  1/8 
-  1/8 
0 

♦  1/8 
-1/8 
0 

♦  1/8 
-  1/8 
0 

♦  !/'l 
0 

♦  1/8 


0  1  0  0. 

1  3.10?  0. 3/18  0. 

1  3.16?  0.3218  0 

1  3.  16?  0.3218  0 

1  3.10?  0.3218  0 

1  3. 162  0.32 18  0 

3  1  3.162  0.3218  0 

3  1  3.  16?  0.3218  0 

'  1  3.162  0.3218  0 

O  1  1.8/08  1. 

3.162  1.86/8  0, 

1 

3.16?  1 .8V? 3  0 

3.16?  1.8923  0 

3.16?  1.8923  0 

3.1416  1. 


s(n  .) 

V.. 

8‘Uc) 

78.608 

100 

0 

0 

1 

1.282 

26.296 

33.709 

27.294 

34.9°  ] 

39 . 944 

51.077 

41.944 

5?. 769 

derivati\e.  The  radius  p(n)  of  the  circle  follows  from  Fig.  19 


As(n)  -  As(n)  for  As  (n)^As  (n- 1 ) 
As(n)  »  is(n-l)  for  A  s  (n)=  As  (n  - 1 ) 


(32) 


The  curvature  x(n)  is  defined  as  the  reciprocal  of  p(n) : 
x(n)  *  1/ p(n)  (S3) 

We  could  draw  circles  according  to  Fig.  19  to  smooth  the 
corners  of  the  polygon  of  Fig.  15,  but  -  apart  from  the  labor 
involved  -  the  approximation  would  become  already  so  good  that 
deviations  would  be  difficult  to  recognize.  Hence,  we  use  a 
different  representation  from  here  on  that  not  only  is  easier 
to  plot  but  will  conspicuously  show  deviations  even  at  the  last 
step  when  the  second  derivative  is  made  continuous. 

Figure  20  shows  on  top  the  sum\Ay(n)  /  2rc  according  to 
column  5  of  Table  7.  The  second  row  shows  Ay(n)/2Tt  according 
to  column  4,  and  the  third  row  shows  Ae(n)  according  to  column 
11.  The  fourth  row  shows  x* (n) = (40 . 58/ 50)x(n)  according  to 
Fqs.(32)  and  (33).  The  star  indicates  that  the  curvature  holds 
for  a  curve  normalized  to  the  length  lOOAx  while  the  length  of 
the  curve  in  Fig.  15  is  2x40. 58Ax  =  81.16Ax. 

Refer  to  Table  7  for  the  computation  of  x*(n).  First, 
\s(n)  is  determined  and  multiplied  by  50/40.58  to  produce  the 
normalized  length  As*(n^.)  listed  in  column  14.  The  normalized 
curvature  x*(n)  follows  then  from  Eqs.(32)  and  (33).  Referring 
now  to  Fig.  19  we  sec  that  the  curvature  calculated  for  the 
point  n  at  the  distance  s*(n)  from  the  arbitrarily  chosen 
beginning  of  the  curve,  also  applies  to  the  interval  with  dis¬ 
tance  +As*(n)/2  from  the  location  s*(n).  We  may  write  for 
this  interval: 

s*  (n)-As*  (n)/2  ^  s*  Si  s*  (n)+As*  (n)/2  (34) 
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Note  that  s*(n)  means  the  normalized  distance  of  the  point  n 
from  the  point  1  measured  along  the  polygon  of  Fig.  15,  while 
s*  means  any  distance  from  the  point  1.  The  difference  is 
that  s*(n)  is  only  defined  for  the  corner  point*  n^. ,  while 
there  is  no  such  restriction  on  s*.  Furthermore,  s *(n)  is  a 
function  of  the  independent  variable  n,  while  s*  i_s  the  new 
independent  variable. 

The  intervals  in  which  the  values  of  x*(n)  apply  are 
shown  by  dashed  lines  in  Fig.  20,  row  4.  For  instance,  the 
value  x*(n)=0.083  in  the  point  s* (nc ) - 1 l . 09  applies  to  the 
range  from  11.09-3.896/2  *  9.14  to  11.09  ♦  3.896/2  -  13.04. 

This  range  is  indicated  in  Fig.  20  by  the  dashed  line  from 
9.14  to  13.04.  The  value  of  this  line  is  given  as  0.083q  for 
the  following  reason.  For  a  closed  curve  the  total  rotation  of 
the  tangent  between  a  point  and  the  same  point  after  advancing 
around  the  curve  equals  2rc,  as  may  be  inferred  from  Fig.  21; 
for  curve;  without  continuous  derivative  one  must  not  only 
start  and  stop  at  the  same  point  but  use  the  same  tangent 
(right  or  left  tangent)  at  the  start  and  the  stop.  For  the 
usual  definition  of  curvature  as  de/ds  we  obtain  the  "total 
curvature"  over  the  length  S  of  a  closed  curve 

/  ds  =  2tt  (35) 

0 

Our  definition  of  curvature  in  F,qs.(32)  and  (33)  must  also 
yield  the  total  curvature  2n  for  a  closed  curve.  We  have 
arbitrarily  assigned  the  normalized  length  100  or  100Ax  to  all 
our  curves,  which  implies  that  we  have  a  value  for  the  curve,  it 
derivatives,  and  its  curvature  at  the  points  s*“0,  l,  ...,  100. 
Hence  we  show  in  Fig.  20,  row  4,  the  value  0.083q  for  the  curva¬ 
ture  at  the  points  s*=10,  11,  12,  13.  Correspondingly,  the  value 
0.052q  is  assigned  to  s*«14,  15,  16;  etc.  The  value  of  q 
follows  from  the  requirement  that  the  total  curvature  must 
equal  2rt  over  its  total  length  S-lOOAx.  From  Fig.  20  -  which 
shows  onlv  half  the  length  S/2  -  follows: 


2(2  x  0.083  ♦  4  x  0.083  ♦  3  x  0.052  ♦  3  x  0.118  +  3  x  0.118  ♦ 

2  x  0.171  ♦  2  x  0.171  +  3  x  0.118  ♦  3  x  0.118  ♦  3  X  0.052 

♦  2  x  0 . 083 ] q  =  2 rt 

q  -  1-02  (36) 

Note  that  the  curvature  in  the  point  0  as  well  as  in  the  point 
100  has  to  be  included,  otherwise  one  would  use  the  right  tan¬ 
gent  at  point  0  but  the  left  tangent  at  the  point  100,  which 
coincides  with  point  0. 

7.  CONTINUOUS  SECOND  DERIVATIVE 

By  rounding  the  corners  in  l-'ig.  15  by  means  of  circles 
according  to  Fig.  19  we  have  made  the  curve  and  its  first  deri¬ 
vative  continuous.  The  curvature  -  or  the  second  derivative  - 
are  not  yet  continuous.  One  may  readily  see  from  Fig.  19  that 
the  points  on  the  circular  section  have  the  curvature  x*  =  l/p*, 
while  the  points  beyond  the  circular  section  have  the  curvature 

x*  =  0. 

In  order  to  make  the  second  derivative  and  the  curvature 
continuous  we  must  develop  an  approximation  that  avoids  the 
discontinuity  of  the  curvature  of  Fig.  19,  but  which  still  is 
based  on  Fig.  19.  The  reason  is  that  the  approximation  of 
Fig.  19  will  detect  circular  contours,  just  as  the  approximation 
by  straight  lines  detected  straight-line  contours.  The  problem 
does  not  have  a  unique  solution.  We  reduce  the  range  of  possible 
solutions  by  further  requiring  that  the  approximation  should  be 
easy  to  implement  by  computer. 

Before  proceeding,  we  must  define  the  normalized  length  of 
the  curve  s*=s*(nc)  in  a  way  that  is  not  changed  every  time  the 
approximation  is  improved,  but  still  identifies  a  point  unambigu¬ 
ously.  In  Figs.  15  and  18  the  length  of  the  contour  was  measured 
along  the  solid  line  connecting  the  corner  points  n  .  The 
approximation  of  Fig.  19  rounded  the  corners  and  thus  shortened 
the  length  of  the  contour,  and  any  further  improved  approxima¬ 
tion  will  change  the  length  of  the  contour  again.  According  to 
Fig.  22  wo  can  use  the  distance  s*  measured  along  the  straight- 
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*  * 

line  approximation  by  connecting  points  Sj,,  ,  ...  on  the 

straight  lines  with  the  centers  P. ,  P,,  ...  of  the  circles  of 

curvature  in  the  points  s^.,  s.,c,  ...  on  the  curved  line.  This 

requires,  of  course,  that  the  centers  P.  ,  P0,  ...  are  known 
*  *  i  * 

before  Sjc>  s-,c>  •••  are  known.  This  requirement  is  satisfied  in 

Fig.  19,  since  the  center  of  the  circle  smoothing  the  break  at 
the  point  n  is  constructed  solely  with  the  help  of  the  straight 
lines  from  n-1  to  n  and  from  n  to  n+1.  Hence,  we  have  the  addi¬ 
tional  requirement  that  the  points  Pj  ,  P^,  ...  in  Fig.  22  must 
be  determined  by  the  straight  lines  from  n-J  to  n  and  from  n  to 
n«- 1 . 

The  approximation  of  Fig.  19  leaves  a  straight  section  with 
curvature  0  at  either  the  longer  one  of  the  lines  from  n-1  to 
n  and  from  n  to  n+1,  or  on  both  lines.  We  look  first  for  a 
method  to  eliminate  these  straight  sections.  According  to  Fig. 

23  we  may  do  so  by  determining  first  a  radius  p^(n)=p(n)  accor¬ 
ding  to  Fig.  19,  and  then  define  a  variable  value  for  p(s)  by 
the  equation 

P(s)  =  P^n)  +  II [ s -  s (n- 1 )  -  As(n-l)/2]  (37) 

where  II  is  a  yet  to  be  determined  constant,  and  s-s(n-l)  - 
As(n-l)/2  is  the  distance  to  the  point  s  measured  from  the  point 
s(n-l)  ♦  \s(n-l)/2.  For  s=s(n+l)  one  obtains  p(  s)  =  P(n  +  1 )  =  t  2  (n) 
and  a  relation  for  II: 


(38) 


If  t he  an;le  6  is  small  compared  with  e(n)  one  obtains  all  the 
relations  necessary  to  determine  Pj (n)  and  P2(n): 

Pj(n)  =  As(n)/2  tan[As(n)/2  -6]  (39) 

P2(m  =  As(n)/2  tan[Ae(n)/2  ♦  p  ] 


P 


*s(n- 1 

TfiT-T 


-  As  fn 
♦  A  s  (n 


tan  Ae(n)/2 


There  are  two  drawbacks  to  this  method.  First,  we  had  to  make 
the  assumption  f><KAr(n)  to  avoid  the  solution  of  a  transcendental 
equation,  which  introduces  either  an  avoidable  inaccuracy  or 
calls  for  a  lengthy  solution  of  a  transcendental  equation  by 
computer.  A  second  drawback  is  even  more  severe.  Let  Ae(n) 
in  Fig.  23  equal  tt  /  2 .  We  will  then  evidently  obtain  pj  (n)  = 
As(n)/2  and  P2(n)  =  As(n-l)/2.  Hence,  the  shorter  interval 
As (n) <  As (n- 1 )  will  get  the  larger  radius  P2(n)  >  Pj(n)  and 
thus  the  smaller  curvature  l/p2(n)  <  1 / P ^  Cn )  in  its  center.  This 
is  just  the  opposite  of  what  we  would  want  for  a  good  approxima¬ 
tion  . 

A  second  method  to  eliminate  the  straight  sections  left  by 
the  approximation  according  to  Fig.  19  is  shown  in  Fig.  24. 

First,  the  curvature  radius  and  the  center  of  the  circle  of  cur¬ 
vature  is  determined  for  every  corner  point  n,  n  +  1,  n  +  2 . 

This  yields  p(n),  p(n+l),  p(n+2) .  and  the  centers  P(n), 

Pln+1),  P(n+2),  ....  of  the  circles.  These  centers  are  connec¬ 
ted  by  the  dashed  lines  with  length  Ad(n) ,  Ad(n+1),  ....  The 
center  of  the  radius  of  curvature  of  an  intermediate  point, 
located,  e.g.,  between  n  and  n+1  and  having  the  distance 
qAs(n)  from  n,  would  then  be  located  on  the  line  from  P(n) 
to  P(n+1)  at  the  distance  qAd(n)  from  P(n);  the  radius  of  cur¬ 
vature  p(n,q)  would  also  change  linearly  with  q  from  p(n) 
to  Pln+1): 

P(n,q)  =  1 1  -q)  p(n)  +  qp(n+l)  (40) 

This  method  would  yield  a  continuously  varying  curvature  and 
thus  a  continuous  second  derivative.  Its  drawback  is  that  the 
radius  of  curvature  pln,q)  could  never  be  larger  than  either 
p(n)  or  p(n+l).  A  look  at  Fig.  15  shows  that  this  means  the 
curvature  anywhere  along  the  straight  line  between  n=l  and  n=10 
could  not  be  less  than  in  the  points  n=l  or  n=10.  This  feature 
makes  this  approach  to  a  continuous  second  derivative  undesirable. 

A  third  extension  of  Fig.  19  is  shown  in  Fig.  25.  It  is 
assumed  that  As(n-l)  is  smaller  than  As(n).  The  verticals  in 


the  center  of  the  line  from  point  n-1  to  point  n,  and  at  the 
distance  As(n-l)/2  from  point  n  on  the  line  from  point  n  to  point 
n  +  1  yield  the  point  P(n)  and  the  radius  p(n) .  The  vertical  in  the 
center  of  the  line  from  point  n  to  point  n  +  l  yields  then  the  point 
P' (n)  and  the  radius  p' (n)  : 

p'(n)  =  pin)  ♦  As'(n)/tan  Ae(n)  (41) 

A  radius  p(n,d)  at  the  point  with  distance  d  from  the  center 
of  the  line  from  point  n-1  to  point  n,  that  varies  continuously 
from  p(n)  to  p'(n),  may  be  defined  in  analogy  to  Eq.(40): 

p(n,d)  =  ( 1  -  q )  p  (n)  +  qp'(n)  (42) 

_  =  _ _  _  d  __ 

q  As(n-l)  A  s'  (nT " 

Substituting  q  and  introducing  the  notation  As(n)  to  include  the 
case  where  As(n-l)  >As(n)  rather  than  As(n-l)  <  As(n)  as  in 
Fig.  25,  one  obtains: 

p(n,d)  =  p(s)  =  P(n)  +  dAs '  (n) /[  As  (.*  /  +  As'(n)]  tan  Ae(n) 

As' (n)  =  [  As (n)  -  As(n-l) |/2  (42) 

As(n)  =  As(n-l),  Sg(n)  =  s(h)  -  As(n)/2  for  As(n-l)  =  As(n) 

=  A s ( n ) ,  s0(n)  =  s(n)  +  As(n)/2  for  As(n-l)  =  As(n) 

d  =  |s  -  Sq  (n)| 

It  is  readily  apparent  that  the  construction  of  Fig.  25  elim¬ 
inates  all  flat  spots  along  a  polygon  as  shown  in  Fig.  15.  The 
radius  p' In)  can  be  arbitrarily  much  larger  than  the  radius  p(n), 
if  As(n)  is  arbitrarily  much  larger  than  As(n-l).  The  curves 
constructed  in  this  way  will  have  a  continuously  varying  curva¬ 
ture,  except  in  the  center  of  the  lines  between  the  points  n-1 
to  n,  n  to  n+1,  etc.  These  points  with  a  discontinuous  second 


derivative  are  a  drawback,  but  in  our  finite  mathematics  one  can 
readily  overcome  this  drawback.  Let  us  explain  the  whole  process 
by  means  of  an  example. 

Table  9  lists  the  values  n,  s*,  Ae  ,  As*(n),  As*(n-1),  As*(n) 
and  p*(n)  =  l/x*(n)  taken  from  Table  7.  Note  that  a  star  charac¬ 
terizes  values  applying  to  the  normalized  length  100  of  the  con¬ 
tour  rather  than  the  actual  length  81.16  of  the  polygon  shown  by 
the  solid  line  in  fig.  15: 

As* (n)  =  (100/81.16)  As(n)  (46) 

The  values  of  s*  in  Table  7  were  only  listed  for  the  corner  points 

n  =1,  10,  13,  ...  while  a  value  of  s*  is  listed  for  every  point 

n  in  Table  9.  These  intermediate  values  are  produced  by  linear 
interpolation.  Let  s*(nc^)  he  the  distance  to  the  corner  point  nc^, 
and  s*(nc2)  the  distance  to  the  next  corner  point  nc2.  The  dis¬ 
tance  s*(n)  to  a  point  n  between  n^  and  nc2  is  defined  by: 


=  s*  (n)  =  s* (n  , )  + 


n  ,  2  n  s  n  0 
cl  c2 


s*(nc2)  -  s*(ncl) 
nc2  *  ncl 


The  value  of  p* (n)  in  Table  9  is  the  reciprocal  of  x*(n)  in 
Table  7.  The  values  of  As'*(n)  and  Sq*  follow  from  F.q.(43).  The 
distances  d  are  divided  into  distances  d^  increasing  with  n  and 
distances  d2  decreasing  with  n: 

dx*  =  s*  -  s0*,  d2*  =  sQ*  -  s*  (45) 

The  radius  of  curvature  p*(n,d)  follows  from  Eq.(43)  and  the  cur¬ 
vature  x*(n,d)  is  its  reciprocal.  We  still  have  to  normalize 
X*(n,d)  in  analogy  to  Fq.(36): 


2q  J1  x*(n,d)  =  5.87q  =  2n 
n^l 

q/ 2k  =  0.170 


Y*/  K  =  0.170  x*(n,d) 
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The  normalized  curvature  m.*/ 2  is  listed  in  column  14  of 
Table  9.  It  is  also  shown  by  the  dots  in  the  lowest  row  of  Fig. 
20. 

In  order  to  obtain  x£/2  not  as  function  of  the  points  n 
but  of  the  normalized  distance  s*  along  the  curve  for  s*  =  0,  1, 

2,  100  -  or  rather  s*  =  0,  Ax,  2Ax . lOOAx  -  we  must 

plot  x£/2rt  for  the  integer  values  of  s*  rather  than  the  integer 
values  of  n  in  Fig.  20.  There  are  many  ways  to  do  so.  The 
simplest  one  is  to  give  r*(s*)  the  same  value  as  x*(n)  has  for 
the  value  of  n  closest  to  s*.  The  values  of  R*/2rt  plotted  for 
s*  =  0,  1,  ...  in  Fig.  20  have  been  chosen  in  this  way.  This 
method  of  determination  of  x*(s*)  from  x*(n)  is  shown  in  more 
detail  in  Fig.  26a. 

A  more  sophisticated  method  that  yields  a  cont inuous  func¬ 
tion  from  the  samples  x*(n)  is  shown  in  Fig.  26b.  The 

samples  for  n  =  1,  2,  ...,  m,  ...  arc  connected  by  the 

dashed  line.  This  yields  -  except  for  a  normalizing  factor  - 
**(s*  )  as  a  continuous  function,  and  we  have  thus  solved  the 
problem  to  derive  from  the  sampled  function  of  Fig.  14a  a  con¬ 
tinuous  function  with  continuous  first  and  second  derivatives. 

The  normalizing  factor  for  the  continuous  function  x*(s*) 
is,  of  course,  zero.  This  causes  no  difficulty,  since  we  will 
never  use  more  than  a  finite  number  of  samples  of  h*(s*).  The 
sample  for  s*  =  j  in  Fig.  26a  follows  readily  from  the  formula 
for  linear  interpolation: 


X*  ( s* ) 


x*(m+l)  -  x  J  (m) 
V,n)  +  s*  (m  +  o  :  s*  (ml 


(48) 


One  must  be  careful  that  another  normalization  is  required  after 
one  makes  the  transition  from  H*  =  x*(n)  to  The 

ellipse  of  Fig.  6  has  72  points  denoted  n,  but  its  normalized 
length  equals  lOOAx,  and  there  are  thus  100  values  x*(s*)  but 
only  36  values  x*(n).  Hence,  we  need  another  normalization 
coef f  ic iei  t 
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(49) 


99 

V1 


>0  £ 


x*(s*)  -  2n 


s"*  0 


The  curvature  with  all  normalizations  is  denoted  x 


N' 


XN  -  Vs*)  =  q0 


(50) 


Note  that  Eqs. (49)  and  (50)  apply  to  the  "step  function"  x*(s*) 
of  Fig.  26a,  to  the  linearly  interpolated  function  x*(s*)  of 
Fig.  2ob,  as  well  as  to  any  function  x*(s*)  derived  from  the 
va  ues  x£(n)  by  a  higher  order  interpolation. 

Fiither  the  step  function  x*(s*)  of  Fig.  26a  or  the  linearly 
interpolated  function  x*(s*)  of  Fig.  26b  can  be  used  to  recog¬ 
nize  the  shape  of  the  contour  in  Fig.  6,  independent  of  shift, 
rotation  or  scaling.  If  one  looks  for  a  contour  of  a  particular 
given  shape,  one  can  cross -corre late  x*(s*)  of  the  observed 
contour  with  x*(s*)  of  the  sample  contours.  This  correlation 
process  provides  an  averaging  over  all  the  observed  values 
Hj(s*)  and  thus  helps  eliminate  the  influence  of  the  square 
field  in  lig.  6.  In  other  cases  one  may  want  to  classify  the 
observed  contour  in  a  more  general  way  than  permitted  by  cross¬ 
correlation  with  sample  functions.  This  will  generally  call  for 
some  averaging  process  to  reduce  further  the  influence  of  the 
square  field  in  Fig.  6.  We  will  show  by  means  of  an  example  how 
to  proceed. 

The  integer  values  of  s*  for  one  half  of  the  observed  ellipse 
are  written  in  the  first  and  sixth  columns  of  Table  10.  The 
values  of  x^/2rt  written  for  integer  values  of  n  but  non-integer 
values  of  s*  in  Table  9  arc  assigned  to  the  integer  values  of 
s*  in  Table  10  according  to  the  staircase  approximation  of 
Fig.  26a.  The  average  values  x*/2n  =  x£(s*)/2tt  arc  obtained 

by  averaging  over  the  11  values  x*(s*'5)/2tt  to  x*(s*+5)/2n  : 


x*(s*)/2rt  =  T|  J  x* (s*  +  i)/2rt 
°  i= -  5  ” 

Summation  >f  the  50  values  of  x*(s*)/2n  for  s*  =  0 


(51) 


4  9  y i e l d  s 
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0.703.  With  the  help  of  Eq. (49)  one  obtains: 

99 

q0  E  X*  (s*)/2n  =  2  X  0.703  q0  «  1  (52) 

s"*0 

qQ  =  0.711 

xn/2tt  =  0.711  x*/2n 

The  values  of  x.,/ 2  are  listed  in  columns  4  and  9  of  Table  10. 

N 

Instead  of  plotting  xn/2tt  we  subtract  first  the  curvature 
of  a  circle  with  length  s*  =  100.  Erom  2TtR  =  100  follows  the 
curvature  x^  =  l/R  =  2rt  / 1 0  0 ,  and  x^/2tt  =  0.01.  This  value 
0.01  is  subtracted  in  columns  5  and  10  of  Table  10.  The  values 
x^/2n-0.01  are  plotted  in  Fig.  27a. 

One  may  readily  recognize  a  typical  feature  of  the  ellipse. 
For  a  short  section  of  the  length  of  the  ellipse  the  curvature 
is  much  larger  than  that  of  a  circle  of  equal  length,  while  for 
a  larger  section  the  curvature  is  smaller  than  that  of  a  circle. 
The  magnitude  of  the  positive  peaks  in  Eig.  27a  thus  exceeds  the 
magnitude  of  the  negative  peaks,  but  there  are  more  negative 
peaks.  The  sum  of  all  positive  and  negative  values  x^/2  -0.01 
is,  of  course,  zero,  since  the  total  curvature  of  any  closed 
contour  equals  2n  according  to  Fig.  21. 

8.  EXAMPLES  OF  CONTOURS 

Let  as  consider  a  few  examples.  First,  we  rotate  the  ellipse 
of  Fig.  C>  into  the  position  shown  in  l’ig.  28.  We  then  produce  a 
table  in  inalogy  to  Table  7,  which  yields  the  corner  points 
n  *  1,  9.  12,  14,  17,  20,  24,  31,  35,  37,  45,  48,  50,  S3,  56, 

60,  67,  71,  as  well  as  the  other  values  corresponding  to  Table  7. 
The  straight-line  approximation  of  the  ellipse  is  shown  in  Fig. 

29.  A  comparison  with  Fig.  15  shows  a  similarity,  but  it  is 
obvious  that  a  much  larger  number  of  corner  points  would  bo  needed 
to  make  Figs.  15  and  29  "approximately  equal". 

Figure  30  shows  VAy(i),  Ay(n) ,  Ae(n),  x*(n)  and  x*/2n  in 
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analogy  to  Fig.  20.  Using  the  staircase  approximation  of  Fig. 

26a,  just  as  it  was  done  in  Fig.  20,  yields  the  normalized  cur¬ 
vature  averaged  over  11  samples  with  the  curvature  of  the  circle 
subtracted  out  as  shown  in  Fig.  27b.  This  Fig.  27b  has  been 
shifted  relative  to  Fig.  27a  so  that  the  maxima  at  n-23  and  n-S9 
in  Fig.  27a  line  up  with  the  maxima  at  n-3l  and  n«6  7  in  Fig. 

27b.  The  location  of  the  maxima  in  terms  of  the  length  s*  of 
the  curve  is  s  *  =  30  and  sr*  *  80  in  Fig.  24a,  and  s^*  *  41 
and  sb*  ■  91  in  Fig.  27b.  The  difference  s^*  -  sR*  ■  11  expresses 
the  angle  of  rotation  vp(b)  between  the  ellipses  of  Figs.  6  and 
28.  This  relationship  between  s^*  -  s,4*  and  vp(b)  is  shown  in 
Fig.  31. 

A  much  better  way  to  obtain  the  angle  v(b)  is  to  cross  - 
correlate  the  functions  of  Figs.  27a  and  b.  The  maximum  of  the 
cross-correlation  function  will  occur  for  a  certain  shift  between 
the  scales  s*,  and  this  shift  is  used  instead  of  the  difference 
sb*  *  s,  *  to  determine  <p(b).  Still  better  is  to  cross-correlate 
the  functions  x*/2tt  in  Figs.  20  and  30. 

Let  us  try  next  a  larger  ellipse  as  shown  in  Fig.  32.  The 
number  of  points  on  its  contour  is  118  rather  than  72  as  in 
Figs,  b  and  28.  The  straight-line  approximation  in  Fig.  33  looks 
much  more  like  an  ellipse  than  the  one  in  Fig.  29.  The  functions 
in  Fig.  34  are  obtained  in  exactly  the  same  way  as  the  ones  in 
Figs.  20  and  30,  except  that  the  linear  interpolation  of  Fig.  2t»b 
rather  than  the  staircase  approximation  of  Fig.  26a  is  used  for 
x*/2n.  1  he  curvature  averaged  over  11  samples  less  the  curva- 

ture  of  the  circle  is  plotted  in  Fig.  27c.  This  plot  is  disap¬ 
pointing  compared  with  the  approximation  of  Fig.  33.  However, 
one  has  to  keep  in  mind  that  Fig.  27c  uses  the  second  derivative 
of  a  continuous  curve  approximating  the  ellipse,  and  the  second 
derivative  of  a  function  is  extremely  sensitive  to  slight  distor¬ 
tions  of  the  function.  We  would  need  an  order  of  magnitude  more 
points  to  make  the  curvature  of  Fig.  27b  look  "smooth",  that  is 
independe  it  of  the  square  field  on  which  the  ellipse  of  Fig.  32 
is  plottet  . 

One  ian  produce  smoother  looking  curvatures  by  averaging 
over  more  than  11  values.  The  average  of  2S  values  of  x£/2tt 
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in  Fig.  20  is  used  -  with  another  minor  modification  -  in  Fig. 

27d.  This  plot  is  much  smoother  than  the  one  of  Fig.  27a  which 
uses  the  average  of  11  values  only.  However,  Fig.  27a  is  better 
since  it  shows  that  there  are  fewer  positive  than  negative 
samples,  and  that  they  have  a  larger  magnitude.  This  characteris¬ 
tic  of  an  ellipse  is  lost  in  Fig.  27d. 

As  the  next  example  we  analyze  the  circular  area  of  Fig.  35. 
The  values  obtained  in  analogy  to  Table  7  are  listed  in  Table  11. 
Due  to  the  symmetry  of  the  circle  we  need  to  list  only  one  quarter 
of  the  points  on  the  contour  -  that  is  76/4  =  19  -  rather  than 
one  half  the  points  as  in  the  case  of  the  ellipse.  One  may  readi¬ 
ly  verify  that  the  circular  line  of  Fig.  36  yields  the  same  values, 
because  the  contour  operators  of  Fig.  7  do  not  distinguish  between 
the  contour  of  an  area  and  a  line.  This  is  so  because  a  line  is 
essentially  a  narrow  area  in  the  mathematics  on  a  square  field. 

One  can,  of  course,  readily  make  a  distinction  by  using  contour 
operators  with  more  than  2x2  points. 

Figure  37  shows  the  straight-line  approximation  defined  by 
columns  6  7,  8  and  11  of  Table  11.  Plots  of^Ay,  Ay,Ae  and 

x*(n)  are  shown  in  Fig.  38  in  analogy  to  Figs.  20,  30  and  34. 

Table  12  carries  the  approximation  further  in  accordance 
with  Table  9.  The  values  of  x*/2ti  are  plotted  at  the  bottom  of 
Fig.  38.  A  linear  approximation  according  to  Fig.  26b  is  used 
to  obtain  x*/2  for  integer  values  of  s*.  These  interpolated 
values  are  listed  in  column  2  of  Table  13,  which  corresponds  to 
Table  10.  The  sliding  average  of  11  samples  of  column  2  is 
listed  in  column  3.  The  sum  of  column  3  yields  0.3317,  which 
means  that  x*/2TC  must  be  multiplied  by  i /  ( 4  x  0.3317)  =  0.7537 
to  obtain  Xn/2tt  in  column  4.  The  normalized  curvature  2rc 
less  the  normalized  curvature  0.01  of  the  "ideal"  circle  is 
listed  in  column  5  and  plotted  in  Fig.  39. 

A  comparison  of  Figs.  27  and  39  shows  a  satisfying  difference 
between  the  ellipses  and  the  circle.  One  should  observe  that  the 
rounding-off  errors  of  the  computations  in  Tables  11  to  13  strongly 
affect  the  plot  of  Fig.  27.  For  instance,  Ac  is  shown  as  0.3218 
for  n=l  but  as  0.3217  for  n=15  in  Table  11.  Hence,  the  rounding 


43 


affects  the  fourth  digit,  and  -  with  two  exceptions  -  the  values 
in  column  5  of  Tabic  13  differ  from  zero  in  the  fourth  digit  only. 

9.  SUMMARY  OF  THEORETICAL  RESULTS 

We  have  shown  that  contour  operators  or  scanning  patterns 
permit  the  detection  of  a  contour  where  the  shade  of  greyness 
of  a  photo  transparency  changes.  This  detection  is  deterministic 
and  can  readily  be  automated. 

We  have  further  shown  that  the  curvature  of  a  contour  or  a 
line  normalized  by  its  length  is  an  invariant  under  shift,  rota¬ 
tion  and  scaling.  A  continuous  normalized  curvature  was  derived 
from  the  finite  samples  taken  with  finite  resolution  Ax  and  Ay, 
or  "on  a  square  field"  for  short.  This  normalized  curvature 
characterizes  the  shape  of  a  contour  or  curve.  The  mathematical 
processes  for  its  derivation  were  chosen  so  that  they  can  be 
readily  automated. 

There  are  two  areas  in  which  the  theory  should  be  carried 
further . 

(a)  The  straight  line  approximations  of  Figs.  33  and  37 
appear  substantially  better  than  the  normalized  curvatures  of 
Figs.  27c  and  39.  The  reason  is  probably  that  the  human  eye 
smooths  the  polygons  of  Figs.  33  and  37  by  looking  at  more  than 
three  adjacent  points,  while  the  theory  developed  so  far  only 
uses  three  (corner)  points.  The  extension  of  the  theory  from 
three  to  more  points  is  straightforward,  and  it  would  probably 
result  in  contour  recognition  and  classification  as  good  as  done 
by  the  human  eye. 

(b)  The  contour  operators  of  Fig.  7  with  2x2  points  each 
are  satisfactory  for  following  a  contour.  They  will  not  distin¬ 
guish  between  an  area  and  a  line  -  as  in  Figs.  35  and  36,  they 
will  not  bridge  a  gap  in  a  line  caused  by  photographic  imperfec¬ 
tion,  and  they  will  not  recognize  the  branching  or  intersecting 
of  lines  -  such  as  road  intersections  and  crossings.  Here  again, 
the  extension  of  the  theory  is  straightforward  to  the  use  of 
contour  operators  with  more  than  2x2  points. 

In  order  to  make  use  of  the  normalized  curvature  for  automated 
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r  all  if  13.  he  'urvutui't  k*  .re  ;‘or  integer  values  ol  s  ,  tie  curvature 
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contour  recognition  one  must  develop  a  library  of  normalized 
curvatures.  One  must  be  able  to  access  them  in  an  organized 
way,  like  one  finds  a  book  in  a  library  by  a  hierarchy  of  more 
and  more  precise  specifications  of  the  subject.  For  instance, 
an  ellipse  would  be  found  by  working  through  the  following 
hierarchy  of  categories: 

1.  straight  line  contour,  yes  or  no 

2.  curvature  of  circle,  yes  or  no 

3.  curvature  periodic  with  half  the  length  of  the  contour, 
yes  or  no 

4.  curvature  equal  to  its  mirror  image  at  four  points,  yes 
or  no 

10.  EQUIPMENT  FOR  CONTOUR  RECOGNITION 
10.1  Block  Diagram 

A  block  diagram  of  equipment  for  contour  recognition  hv 
means  of  contour  operators  is  shown  in  Fig.  40.  The  heart  of  the 
equipment  is  a  plasma  display  tube  with  a  resolution  of  512  x  512 
points.  Such  tubes  are  commercially  available  (Petty  and  Slottow 
1971).  The  tube  can  produce  patterns  like  the  Walsh  or  block 
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pulse  patterns  of  Fig.  8  anywhere  in  its  field  of  512“  points. 

Such  a  pattern  with  8x8  points  -  rather  than  2x2  points  as  in 
Fig.  8  -  is  shown  in  Fig.  40.  The  equipment  is  designed  so  that 
any  pattern  with  n  x  m  points,  where  n,  m  =  1,  2,  4,  8,  32,  64, 
can  be  produced.  This  makes  the  equipment  good  enough  for  the 
generalization  mentioned  in  Section  9  under  point  b. 

The  plasma  tube  display  requires  512  power  drivers  for  the 
variable  x  (PPRx)  and  512  more  for  the  variable  y  (PPRy) .  The 
selector  switches  SSWx  and  SSWy  determine  which  power  drivers 
are  activated.  The  format  switches  FSWx  and  FSWy  determine  the 
format  of  the  scanning  pattern  in  Fig.  6  and  the  contour  operators 
in  Fig.  7.  As  pointed  out, the  equipment  is  designed  for  the  more 
general  format  n  x  m  rather  than  2x2.  The  binary  function 
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generators  BFGx  and  Bl  Gy  produce  functions  according  to  Fig.  8 
hut  with  up  to  64  x  64  point  elements  instead  of  the  2x2 
shown.  A  display  drive  control  DDC  produces  the  commands  for 
the  various  circuits. 

The  basic  idea  for  the  use  of  the  equipment  of  Fig.  40  is  as 
follows:  A  transparency  of  an  aerial  photograph  is  placed  over 

the  plasma  tube  display.  The  display  produces  a  set  of  func¬ 
tions  according  to  Fig.  8  and  scans  step  by  step  in  the  x-  and 
y- direct  ions .  The  light  passing  through  the  transparency  is 
detected  by  a  photo  tube,  the  output  voltage  of  which  represents 
the  coefficients  b(i,k)  of  Eq.(19),  or  similar  coefficients  from 
which  the  Walsh  coefficients  a(i,k)  can  be  derived  as  in  F.q.(20). 

By  digitizing  the  output  voltage  of  the  photo  tube  one  may  thus 
produce  the  Walsh  coefficients  and  use  them  according  to  Table  2 
to  determine  the  observed  contour  operator.  The  observed  contour 
operator  then  determines  according  to  Fig.  7  in  which  direction 
the  scanning  pattern  on  the  plasma  tube  display  should  be  advanced. 
This  yields  the  values  of  and  At)  in  Table  1  or  Table  6,  from 
which  everything  else  is  derived  (See,  e.g.,  Tables  7,  8  and  11). 
Hence,  in  addition  to  the  blocks  shown  in  Fig.  40  we  need  a  photo 
detector  with  digitizer,  and  a  computer. 

A  refinement  of  Fig.  40  is  not  to  place  the  transparency 
directly  on  the  plasma  tube  display,  but  to  project  an  image  of 
the  plasma  tube  display  on  the  transparency  by  means  of  optical 
lenses.  There  arc  two  reasons  for  it.  First,  the  glow  discharge 
inside  the  plasma  tube  is  about  two  millimeters  distant  from  the 
transparency,  which  reduces  the  resolution.  A  lens  system  can 
project  the  glow  on  the  transparency  so  that  it  is  in  focus  there 
and  yields  the  best  resolution.  Second,  the  use  of  the  lens 
system  permits  a  reduction  of  the  projected  image  of  the  plasma 
tube.  The  resolution  of  these  tubes  is  about  2  (glowing)  lines 
per  millimeter.  A  reduction  of  10:1  increases  the  resolution  to 
20  lines  per  millimeter,  which  is  closer  to  the  resolution  of 
photo  transparencies.  A  typical  aerial  photograph  would  require 
about  10  x  10  such  reduced  images  of  the  plasma  tube  to  be 
covered,  which  implies  that  the  transparency  must  be  moved 
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mechanically  10  steps  in  the  x-direction  and  10  steps  in  the  y- 

direction  to  be  fully  analyzed.  There  is  currently  no  way 

? 

around  this.  The  available  plasma  displays  resolve  512“  ■ 

262,144  points,  and  the  mechanical  movement  increases  this  num¬ 
ber  to  512“  x  10“  -  262,144,000  points.  This  is  the  order  of 
magnitude  needed  to  analyze  aerial  photographs.  If  and  when 
electronic  displays  with  that  many  resolved  points  become  avail¬ 
able,  we  can  do  away  with  the  mechanical  movement.  The  equip¬ 
ment  described  below  does  not  have  the  mechanical  stepping 
device  required  to  make  10  steps  in  the  x-direction  and  10  steps 
in  the  y-Jirection,  but  it  is  designed  so  that  this  improvement 
can  be  added  readily. 

In  t  tie  following  sections  we  describe  in  more  detail  15 
sheets  of  blueprints  that  contain  the  mechanical  design,  elec¬ 
tronic  circuits  and  timing  diagrams.  Due  to  their  size  and 
specialized  interest  the  blueprints  are  not  supplied  with  all 
copies  ot  this  report  but  are  added  to  a  few  selected  copies 
only.  Adlitional  sets  of  the  blueprints  may  be  requested  from 
the  author. 

10.2  Meehan i ca 1  Dos i gn 

The  mechanical  design  is  based  on  two  enclosed  telephone 
cabinets  is  shown  in  the  blueprint  Plasma  Image  Analyzer,  Mech¬ 
anical.  'n  the  left  are  the  electronic  circuits  that  drive  the 
plasma  tu >e  display.  On  top  of  them  is  the  plasma  display  itself. 
The  light  from  the  display  is  deflected  by  Mirror  1  into  the 
cabinet  on.  the  right ,  and  by  Mirror  2  down  through  the  lens  assem¬ 
bly  to  the  fresnel  lens.  The  photo  transparency  is  to  be  placed 
on  top  of  the  fresnel  lens,  and  the  image  of  the  plasma  tube 
display  i.  focused  onto  it.  In  its  current  design,  the  image 
of  the  plasma  display  is  reduced  only  slightly  to  fit  the  stan¬ 
dard  format  of  aerial  photographs.  To  achieve  a  reduction  10:1 
one  must  change  the  lens  assembly  and  move  it  closer  to  the 
fresnel  lens.  The  design  makes  this  readily  possible.  further¬ 
more,  a  mechanical  scan  can  easily  be  added  by  stepping  the 
phototran  iparency  ten  steps  in  the  direction  vertical  to  the 
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paper  plane  (y-d i  roc  t  ion )  and  mov  ing  tlie  projected  image  by 
means  of  mirrors  ten  steps  from  left  to  right  in  the  paper  plane 
(x -d ireet ion) . 

The  fresnel  lens  collects  most  of  the  light  passing  through 
the  transparency  and  throws  it  on  a  photo  tube.  The  output 
voltage  of  this  tube,  after  some  processing  and  digitizing,  is 
fed  to  a  computer  which  produces  the  command  what  to  do  next  for 
the  Pi  splay  Drive  Control  1>I)C  in  the  left  telephone  cabinet. 

The  blueprints  Mirror  Mount  1,  Mirror  Mount  2  and  Fresnel 
Lens  Mount  show  details  fo  the  design  of  these  items.  The  adjust¬ 
ment  screws  with  springs  in  the  mirror  mounts  permit  the  fine 
adjustment  of  the  optical  axis. 

10.3  Pr inc ip  1 e  of  electronic  Design 

The  block  diagram  shown  in  the  blueprint  Plasma  Display 
Image  Analyzer,  lilectronic  begins  with  the  computer  shown  in 
the  center.  Through  the  Computer  Interface  Circuit  C1C  one 
accesses  the  Display  Drive  Control  DDC  and  the  Fototube  and 
Sensor  Circuit  SHN.  The  data  flow  is  indicated  by  arrows.  The 
circuit  design  is  based  on  computer  architecture  principles.  A 
number  of  lines  such  as  Data  Bus  and  Address  Bus  are  connected 
to  all  circuits,  and  the  address  sent  on  the  Address  Bus  deter¬ 
mines  which  circuit  will  receive  the  data  on  the  Data  Bus.  Table 
14  lists  the  common  lines.  Table  15  explains  the  functions  of 
CON'l  ;OL  1  and  CONTROL  2  in  more  detail.  Table  16  lists  the 
address  odes  of  the  addressable  circuits. 

The  general  operation  is  as  follows.  The  controlling  com¬ 
puter  sends  an  address  on  the  ADPRF.SS  BUS.  CONTROL  2  enables 
all  circuits  to  receive  this  address,  but  only  the  one  with  the 
correct  address  will  be  activated.  For  instance,  the  address 
1111  will  activate  the  format  switch  control  FSC  for  x.  If  the 
computer  sends  now  a  number,  c.g.,  4=1000,  on  the  DATA  BUS,  the 
format  switch  control  FSC  will  close  the  first  four  format 
switches  I SW.  Cenerally,  each  operation  requires  two  numbers 
from  the  control  computer,  the  address  and  the  data.  IVe  will  see 
later  that  there  is  one  exception:  The  gas  discharge  in  the 


Table  14.  Common  lines  for  all  circuits 
of  the  Plasma  linage  Analyzer 


Name 

Function  | 

1 

CLOCK  1 

main  clock  pulses 

7 

4m 

CLOCK  2 

main  clock  pulses  shifted 
half  a  clock  period 

by 

3 

CONTROL  1 

enables  data  read- in 

4 

CONTROL  2 

enables  address  read  in 

5 

ADDRESS  HUS 

transmits  address 

6 

DATA  BUS 

transmits  data 

7 

1N1IIB1  r 

inhibits  data  and  address 
transfer  while  an  instruct 
is  operated  on 

i  on 

Table  15.  Function  of  CONTROL  1  and 
CONTROL  2  signals. 


CONTROL  1 

Enables  data  read-in.  HIGH  during  data  read-in,  LOW 
otherwise . 

During  HIGH  data  can  be  read  serially  into  the  acti¬ 
vated  shift  register.  The  transition  from  HTGH  to  LOW 
transfers  the  data  paralled  from  these  registers  to 
buffer  storages. 


CONTROL  2 

Enables  address  read-in.  HIGH  during  address  read- 
in,  LOW  otherwise. 

During  HIGH  the  address  code  is  read  into  the  shift 
registers  of  all  address  decoders.  The  pulses  from 
CLOCK  1  are  inhibited  at  all  other  times.  Hence, 
the  last  address  code  transmitted  is  stored  in  all 
shift  registers  of  the  address  decoders  until  the 
next  address  code  arrives. 


Table  16.  Assignment  of  Addresses 


address  number 

addressed  circuit 

0 

0000 

1 

0001 

L 

0010 

3 

0011 

4 

0100 

5 

0101 

6 

0110 

Rene ra 1  b  i  na r y  func t i on  cir¬ 
cuit  for  y 

7 

0111 

general  binary  function  cir¬ 
cuit  for  x 

8 

1000 

voltage  control  for  photo 
mu ltiplier 

9 

1001 

gas  discharge  control 

10 

1010 

display  drive  control  for  y 

l  1 

1011 

display  drive  control  for  x 

1  2 

1100 

selector  switch  control  for  v 

13 

I  101 

selector  switch  control  for  x 

14 

1110 

format  switch  control  for  v 

1  5 

1 1 1 1 

format  switch  control  for  x 

L 
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display  is  switched  on  by  the  transmission  of  an  address  alone. 

The  circuit  shown  in  the  blueprint  Plasma  Image  Analyzer, 
Flectronic  contains  the  same  circuits  for  the  variables  x  and  y. 
Table  lb  shows  the  address  numbers  assigned  to  the  various  cir¬ 
cuits.  Only  the  addresses  6  to  15  are  currently  used,  the  others 
are  a  reserve  for  a  future  expansion  of  the  equipment. 

1'he  Selector  Switches  SSK  and  the  Power  Drivers  PDR  are 
divided  into  8  parts  SSh'l  -  SSW8  and  PDR1  -  PDR8 .  The  reason  is 
that  312  output  terminals  are  required  for  these  circuits  to  drive 
the  512  lines  or  columns  of  the  plasma  tube  display,  and  it  is 
very  difficult  to  accomodate  that  many  terminals  on  fewer  than 
eight  printed  circuit  cards. 

The  various  circuits  shown  in  the  blueprint  Plasma  Image 
Analyzer,  Flectronic  are  discussed  in  more  detail  in  the  follow¬ 
ing  sections. 

10.4  Display  Drive  Control  PDC 

Refer  to  the  blueprint  Display  Drive  Control  DDC.  It  con¬ 
tains  three  sections.  One  produces  the  function  wal(k,x),  the 
second  the  function  wal(m,y),  and  the  third  fires  the  plasma 
display.  All  three  sections  have  to  be  addressed  by  a  control 
computer  via  the  ADDRFSS  HUS  and  CONTROL  2. 

Let  the  address  code  1011  be  transmitted  on  the  ADDRFSS  BUS. 
The  aditress  decoder  11  for  the  function  wal(k,x)  will  be  acti¬ 
vated,  which  means  that  a  HIGH  voltage  is  obtained  at  terminal  \\' 
of  MUL3  if  a  HIGH  voltage  is  applied  to  CONTROL  1.  This  enables 
the  shift  register  SRF5,  and  the  number  on  the  DATA  BUS  is  read 
into  SRF5.  The  number  is  transferred  parallel  to  the  D-type  flip- 
flop  storage  I DH2  at  the  end  of  the  enable  signal  from  MUL3. 

This  number  represents  the  normalized  scquency  k  of  wal(k,.\). 

It  is  used  to  produce  the  spatial  function  wal(k,x)  in  the  way  to 
be  discussed  in  the  following  Section  10.5. 

The  address  code  1010  on  the  ADDRFSS  BUS  will  activate  t  lie 
address  decoder  10  and  produce  the  spatial  function  wal(m,x)  in 
complete  analogy  to  the  generation  of  wal(k,x). 

The  address  code  1000  on  the  APURFSS  BUS  will  activate  the 


address  decoder  9.  Every  time  this  address  is  received,  the  plasma 
display  is  fired  twice;  during  the  first  discharge  half  of  the 
wires  of  the  plasma  display  selected  by  the  Format  Switches  will 
be  connected  to  a  positive  voltage,  the  other  half  to  a  negative 
voltage,  while  the  polarity  is  reversed  during  the  second  dis¬ 
charge.  Four  output  signals  have  to  be  produced  by  the  Timer  for 
Power  Drivers  (TPD)  to  generate  this  double  discharge.  For  an 
explanation  of  this  circuit  refer  to  the  blueprint  Time  Diagram 
for  TPD,  which  shows  the  voltages  at  the  points  a,  b,  ...,  p  as 
function  of  time.  The  enable  signal  at  point  a  is  a  positive 
step  function.  It  is  differentiated  to  produce  a  positive  trigger 
pulse  at  point  b,  and  the  same  pulse  with  reversed  amplitude  at 
point  c.  A  timing  circuit  TIM1  produces  a  positive  pulse  of 
about  20  ps  duration  at  point  d;  the  duration  of  this  pulse  can 
be  adjusted  as  indicated  at  the  falling  edge  of  the  pulse  in 
line  d  of  the  Time  Diagram  for  TPD.  A  modulo  2  counter  FFD  pro¬ 
duces  al ternat ingly  a  HIGH  and  a  LOW  voltage  at  point  e  every  time 
a  positive  edge  is  produced  by  TIM1.  The  gate  AND  2.1  passes  the 
voltage  from  point  d  whenever  point  e  is  HIGH.  Note  the  short 
pulse  (glitch)  in  line  f  of  the  Time  Diagram  for  TPD;  it  must  be 
eliminated  for  a  clean  operation  of  the  circuit. 

The  timer  TIM2  produces  a  pulse  of  2  ps  duration  at  point  g 
for  every  negative  edge  at  point  f.  This  pulse  is  inverted  at 
point  h,  and  differentiated  at  point  i;  the  positive  pulse  at 
point  i  bi'comes  a  negative  pulse  at  point  c,  while  the  negative 
pulse  at  ooint  i  is  suppressed.  The  action  of  the  negative  pulse 
at  the  points  c  to  i  is  shown  in  the  blueprint  Time  Diagram  for 
TPD  by  dashed  lines.  The  main  feature  is  that  the  modulo  2  coun¬ 
ter  1 FD  turns  off  the  gate  AND  2.1  and  thus  prevents  a  second 
cycle  of  the  feedback  loop  to  reach  point  i.  Hence,  the  enable 
voltage  at  point  a  produces  the  voltages  at  the  point  j  and  k 
shown  in  the  Time  Diagram  for  TPD  in  the  interval  0  <  t  <  T,  and 
no  more. 

I  he  gate  NOR2  produces  negative  trigger  pulses  at  point  l, 
whic  are  transformed  into  rectangular  pulses  of  1  |is  duration 
at  p  int  n  by  the  timer  TIM3;  the  inverted  pulses  are  produced 


at  point  n.  The  signal  with  a  glitch  is  produced  at  point  o, 
the  glitch  is  suppressed  by  an  RC- filter,  and  clean  signals  are 
derived  at  the  points  p  and  p.  These  two  signals  will  be  needed 
to  turn  on  the  power  drivers  of  the  plasma  display.  The  signals 
e  and  e  are  used  to  produce  first  the  Walsh  functions  -wal(k,x), 
♦wal(m,y)  and  then  their  complements  +wal(k,x),  -wal(m,y).  Note 
that  the  transients  of  e  and  e  occur  when  the  voltage  at  point  p 
is  high,  and  this  prevents  damaging  surge  currents  during 
t  runs  ient  s . 


10.5  Binary  function  Generato r  BIG 

The  Binary  Function  Generator  BIG  comes  in  two  versions. 

The  first  accepts  any  string  of  binary  numbers  from  the  computer 
and  stores  them  as  function  of  either  x  or  y.  This  permits  the 
generation  of  any  binary  function  but  the  process  of  loading  is 
time-consuming.  A  second  version  produces  only  the  Walsh  func¬ 
tions,  but  it  is  hard  wired  and  thus  very  fast. 

For  a  discussion  of  the  first  version  refer  to  the  blueprint 


General  Binary  Function  Circuit.  Upon  reception  of  the  address 
code  6  for  y  and  7  for  x  the  digits  on  the  DATA  BUS  will  be  stored 
in  the  shift  register  SRfc'l  to  SRE8 .  A  total  of  64  digits  can 
be  stored,  but  only  2  digits  are  needed  to  produce  the  four  block 
functions  of  Fig.  8.  For  such  a  small  number  of  digits  the  cir¬ 
cuit  works  about  as  fast  as  the  hard  wired  Walsh  function  circuit. 


The  inverters  1NV3  and  INV4  delay  the  clock  pulse  to  make  sure 
the  voltage  at  terminal  W  has  reached  the  HIGH  state  when  the 
clock  pulse  reaches  the  shift  registers  SRE1  to  SRF8 .  At  the 
end  of  the  enable  pulse  the  digits  in  the  shift  registers  are 
transferred  in  parallel  to  the  D-tvpe  flip-flops  FPH1  -  FDHll. 

The  specialized  circuit  for  Walsh  functions  only  is  shown 
in  the  blueprint  Binary  Function  Generator  for  Walsh  Functions 
BIG  and  Format  Switches  FSW.  This  circuit  lias  been  discussed 


in  the  literature  (llarmuth  1977  ,  p.150).  The  voltages  at  the 
input  terminals  xl  to  xo  are  either  HIGH  (1)  or  LOW  (10).  The 
states  of  the  output  terminals  corresponding  to  the  2^  combina¬ 
tions  of  voltages  at  the  input  terminals  arc  shown  in  the 


blueprint  Time  Diagram  for  Binary  Function  Generator  for  Walsh 
Functions. 

The  format  switches  FS  are  field  effect  transistor  switches 
with  a  very  high  impedance  in  the  open  state.  Hence,  the  output 
terminals  l  of  the  format  switches  FS  have  three  states:  HIGH 
(+5V) ,  LOW  (Ov)  and  FLOATING  (no  specified  voltage  but  essentially 
no  current).  The  purpose  of  the  format  switches  is  to  select  the 
format  n  x  m  of  the  functions  used.  All  formats  with  n  =  1,  2,  4, 
8,  16,  32,  64  and  m  =  1,  2,  4,  8,  16,  32,  64  are  possible;  n 
and  m  can  be  selected  independently. 

10.6  Format  and  Selector  Switch  Control 

Refer  to  the  blueprint  Selector  Switch  Control  SSC  and 
Format  Switch  Control  FSC.  The  circuit  diagram  for  the  Format 
Switch  Control  starts  out  with  an  address  decoder  set  for  the 
code  15  =  1111  for  the  variable  x  and  14  =  1110  for  the  variable 
y.  The  correct  address  will  enable  the  read-in  from  the  DATA 
BUS  to  the  shift  register  SRE5  when  the  line  CONTROL  1  is  HIGH. 

A  six  digit  number  is  read  in,  and  transferred  in  parallel  to 
the  D-type  flip-flop  storage  FDH  at  the  end  of  the  enabling 
signal.  The  OR-gates  make  some  or  all  of  the  output  terminals 
n  HIGH  according  to  the  listing  of  Table  17.  The  voltages 
H  =  +5V  and  L  =  0  V  open  or  close  the  format  switches  in  the  blue¬ 
print  Binary  Function  Generator  for  Walsh  Functions  BFG  and  Format 
Switches  FSW. 

The  Selector  Switch  Control  SSC  works  very  similarly.  The 
address  codes  13  =  1101  for  x  and  12  =  1101  for  y  permit  the  read- 
in  from  the  DATA  BUS.  The  data  is  transferred  in  parallel  at 
the  end  of  the  enabling  pulse  to  the  storages  FFD1  and  FFD3.  A 
total  of  9  binary  digits  is  stored,  yielding  29  =  512  different 
output  states.  Both  the  outputs  Q  and  their  complements  Q  are 
used.  The  2x9  voltages  are  fed  to  the  Selector  Switches. 

Each  terminal  Q  or  Q  can  be  loaded  with  16  mA,  but  only  512  x  10 
pA  =  5.12  mA  arc  required  to  drive  the  switches.  Hence,  no 
overloading  occurs  despite  the  512  s i ngle -pole- single  -  throw 
switches  driven  by  each  output  terminal. 

-  58  - 


I 


.  nl'lt'  IV-  1 1 1 J  it  data  on  DATA  HUS  and  states  or  the  output  tormina 
i  of  the  Format.  Switch  Control  FSO.  The  first  digit  of  the  input 
hat  a  arrives  first,  and  is  stored  in  Q^. 
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10.7  Selector  Switches  SSW 

Refer  to  the  blueprint  Principle  of  Selector  Switches.  It 
shows  eight  input  terminals  l  -  1  ...  8,  six  sets  of  64  switches 
each,  and  2^  =  64  output  terminals  i  =  1  ...  64.  All  switches 
are  in  the  position  0  (right);  when  they  are  flipped  to  the  left 
they  are  in  the  position  1.  Let  the  eight  voltages  a,  b,  ...  h 
be  applied  to  the  input  terminals  l.  The  voltages  at  the  output 
terminals  for  various  positions  of  the  six  sets  of  switches  SSl,i 
to  SS6,i  with  i  =  1  ...  64  are  shown  in  Table  18.  The  output 
terminals  that  have  no  voltage  assigned  are  FLOATING;  they  will 
permit  practically  no  current  to  flow. 

The  implementation  of  this  principle  for  up  to  64  input 
voltages  at  the  terminals  l  -  1  ...  64,  nine  sets  of  switches 

9 

SSI ,  i  to  SS9,i  and  2  =  512  output  terminals  i  =  1  ...  512  is 

shown  in  the  blueprint  Selector  Switch  Card  SSC.  Actually  only 
64  output  terminals  i  are  shown  since  one  needs  8  such  cards 
that  are  connected  as  shown  in  the  blueprint  Plasma  Image  Analy¬ 
zer,  Electronic. 

10.8  Power  Drivers  PDR 

Each  selector  switch  card  is  connected  to  a  power  driver 
card.  The  circuit  diagram  of  these  cards  is  shown  in  the  blue¬ 
print  Pow:;r  Driver  Card  PDR.  Each  card  contains  64  power  drivers 
PD,  vhose  circuit  is  shown  in  the  upper  left  corner. 

The  >ower  drivers  produce  the  three  output  states  0,  +125  V 
(  +  12<).  aid  floating  [F]  in  response  to  the  three  input  states 
0,  +>  V  (-5),  and  floating.  The  four  transistors  T3  -  T6  must 
be  high  voltage  types,  such  as  2N5059  or  A5T5059.  The  rather 
larg  ■  num >er  of  four  transistors  is  chosen  to  save  power,  which 
is  m  indat  >ry  when  working  with  the  high  voltages  required  by  a 
plasma  di  .play. 

The  ,.ates  AND  and  OR  prevent  overloading  of  the  transistors 
T5  and  T6  during  transients.  For  an  explanation,  let  T6  be 
conducting  and  T5  by  nonconducting  (input  voltage  0,  output 
voltage  0).  If  the  input  voltage  is  now  changed  from  0  to  +5  V, 
the  trans  stor  T5  becomes  conducting  and  the  transistor  T6 
nonconducting.  During  the  transient,  T5  as  well  as  T6  may  be 
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conducting,  and  the  high  voltage  of  12S  V  could  quickly  cause 
thermal  destruction.  The  gates  OR  and  AND  make  both  15  and  T6 
nonconducting  during  any  transient,  and  only  after  a  certain 
guard  time  will  one  of  the  two  transistors  be  made  conducting. 

The  AND-gate  is  controlled  by  the  voltage  on  the  line  p  of  the 
Timer  for  Power  Drivers  (blueprints  Display  Drive  Control  DDC  and 
Time  Diagram  for  TPD) ,  while  the  OR-gate  is  controlled  by  the 
voltage  on  line  p. 

10.9  Pototube  and  Sensor  Circuit  SEN 

Refer  to  the  blueprint  Pototube  and  Sensor  Circuit  SEN.  The 
current  of  the  fototube  is  converted  to  a  voltage  and  amplified 
to  the  1  V  level.  Via  the  sample-and-hold  circuit  SUM  the  voltage 
is  fed  to  the  analog-digital-converter  ADC,  which  feeds  the 
digits  to  the  control  computer. 

Two  unusual  circuits  are  required.  A  programmable,  protected 
power  supply  adjusts  the  operating  voltage  of  the  fototube  under 
computer  control  to  the  brightness  of  the  plasma  tube  display. 

The  reason  for  this  is  that  the  format  of  the  displayed  two- 
dimensional  function  may  vary  from  mxn=lxl  =  ltomxn  = 

64  x  64  =  4112.  The  light  emitted  by  the  plasma  display  may  thus 
vary  by  4112:1,  and  the  variation  of  the  light  received  by  the 
fototube  is  just  as  large.  One  can  adjust  the  amplification  of  a 
fotomul t ipl ier  tube  by  changing  its  voltage.  This  is  the  purpose 
of  the  programmable  power  supply. 

The  address  8  =  1000  will  enable  the  read- in  of  data  from 
the  DATA  BUS  to  the  shift  register  SRE2 .  At  the  end  of  the 
enabling  pulse  the  data  is  transferred  in  parallel  to  the 
stor  iges  PFD1  and  I  PD2.  The  d ig ital - to -analog  converter  DAC 
transforms  the  digits  into  a  voltage  which  controls  the  output 
voltage  of  the  power  supply  PSM. 

A  fototube  is  easily  overloaded  by  accidental  illumination 
from  other  sources  than  the  plasma  tube  or  by  incorrect  setting 
of  the  operating  voltage  delivered  from  the  programmable  power 
supply.  A  protection  circuit  is  used  that  sets  the  operating 
voltage  to  zero  if  the  fototube  is  saturated  for  about  one  second. 
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To  do  so  the  output  voltage  ot'  0PA4  is  integrated  bv  0PA2.  The 
integrator  is  reset  to  zero  once  a  second  via  the  timer  T1M4. 

If  during  such  an  one  second  interval  the  integrated  voltage 
exceeds  about  10  V,  a  flip-flop  1;JK  is  triggered  which  switches 
the  reference  voltage  of  the  d ig i ta 1 -analog-converter  via  the 
switch  ASW1  from  -4.7  V  to  0.  At  the  same  time,  a  glow  lamp 
indicator  "OUT"  is  turned  on.  The  voltage  for  the  t'ototube  must 
he  reset  manually  by  means  of  t  lie  switch  RESF.T  "OUT". 

11.  CONCLUSION  ANU  RECOMMENDATION 

The  contour  recognition  and  classification  based  on  angles 
and  curvature  probably  simulates  very  closely  the  processes  used 
bv  the  human  eye  and  brain.  It  is  a  much  more  sensitive  method 
than  classifications  using  correlation  and  mean  square  errors, 
which  are  based  on  integration,  since  it  uses  the  second  deri¬ 
vative  of  the  contour  which  is  much  more  sensitive  to  small 
changes  than  an  integral  over  the  area  surrounded  by  the  contour. 

The  equipment  needed  to  implement  the  theory  is  well  within 
the  state  of  the  art.  No  particular  precision,  as  in  Fourier 
optics  processing,  is  needed,  and  the  equipment  lends  itself 
readily  to  computer  control  in  a  feedback  loop. 

The  following  recommendations  are  made  for  the  next  steps  of 
this  program: 

(a)  Equipment  according  to  the  blueprints  should  be  developed. 
The  cost  would  be  in  the  order  of  $100,000  ,  if  the  development 
is  done  by  someone  with  previous  experience  with  pjasma  tube  dis¬ 
plays  and  a  scientific  interest  in  the  project,  while  about 
twice  this  amount  will  be  required  by  a  reliable  organization  with 
purely  financial  interest  and  no  previous  experience  with  plasma 
tubes.  Before  proceeding  with  the  development  a  final  de tormina 
t ion  should  be  made  how  versatile  and  complex  one  wants  the  equip¬ 
ment  to  be,  since  t be  development  of  the  theory  of  contour  recog¬ 
nition  based  on  angles  and  curvature  permits  considerable  simplifi¬ 
cations  compared  with  the  equipment  originally  called  for  in  this 
contract.  After  this  possible  final  modification  the  design 


should  be  frozen,  since  there  arc  as  many  variations  of  the  design 
as  there  are  engineers  and  the  perfect  equipment  is  never  built. 

(b)  There  are  three  directions  in  which  the  theory  should 
he  expanded,  but  this  advancement  of  theory  is  only  justified  if 
equipment  is  developed.  The  first  direction  is  to  work  out  a 
continuous  second  derivative  of  a  contour  which  is  based  on  more 
than  three  adjacent  points;  this  appears  to  be  at  this  time  the 
primary  shortcoming  of  the  theory  compared  with  a  human  observer. 
The  second  direction  is  to  work  out  an  intelligently  organized 
"library"  of  typical  contours  encountered  in  aerial  photogtaphy, 
which  is  needed  to  be  able  to  classify  the  normalized  curvature 
of  an  observed  contour.  The  third  direction  is  to  advance  from 
scanning  patterns  and  contour  operators  with  2x2  points  to  more 
points,  in  order  to  solve  the  problem  of  road  crossings  and 
related  problems. 


12.  COMMENTS  ON  TF.CIINOl.OC.Y 


An  aerial  photograph  resolves  about  10(1(1(1“  points.  In  order 
to  implement  the  concepts  in  this  report  one  needs  an  electro¬ 
optic  converter  that  potentially  permits  this  resolution,  can 
operate  with  a  computer  in  a  feedback  loop,  and  is  sufficiently 
fast. 

The  cathode  ray  tube  is  the  most  widely  used  electro-optic 

converter.  It  is  readily  available  for  the  resolution  of  about 
"> 

1000“  points.  However,  the  existing  diplays  have  been  developed 
for  applications  with  very  different  requirements  and  are  thus 
not  well  suited  for  image  processing.  For  instance,  even  though 
the  writing  speed  of  an  electron  beam  is  in  the  order  of  centi¬ 
meter  per  nanosecond,  the  decay  time  of  the  phosphor  is  purpo¬ 
sely  made  very  long,  typically  in  the  order  of  milliseconds  and 
longt  r.  Tiiis  makes  such  a  converter  as  slow  as  a  liquid  crystal 
display  i  i  this  application.  A  second  inherent  problem  is  the 
acctii  icv  <  f  the  absolute  location  of  the  electron  beam.  A  cathode 
ray  tube  1  av  write  1000  lines  on  a  TV  monitor,  but  these  lines 
may  be  sh  fted  up  or  down  In-  a  few  line  widths  without  detrimen¬ 
tal  effect,  as  long  as  all  the  lines  are  shifted  equally.  This  is 
not  tolerable  for  image  analysis,  since  one  would  not  know  exactly 
which  point  of  an  image  is  hoeing  analyzed. 

The  problem  of  absolute  location  of  the  cathode  ray  tube  is 
overcome  in  the  crossbar  systems,  the  best  known  of  which  are 
liquid  crystal  displays  and  plasma  displays.  The  crossbars,  con¬ 
sisting  of  wires  or  conducting  strips,  are  mechanically  rigid 
and  cannot  shift  due  to  voltage  fluctuations.  The  main  problem 
of  crossbar  systems  is  resolution. 

The  largest  liquid  crystal  display  known  to  the  author  re¬ 
solve-  32'  points,  and  no  work  has  been  reported  recently  about 
an  increase  of  this  resolution.  Furthermore,  the  liquid  crystal 
display  i?  slow,  requiring  upwards  of  1  ms  to  advance  from  one 


pattern  to  the  next.  Among  its  good  features  are  the  excellent 
definition  of  the  patterns,  low  driving  voltage,  and  next  to 
no  driving  power. 

2 

Plasma  displays  are  available  with  a  resolution  of  512 

points.  They  can  switch  from  one  pattern  to  the  next  in  about 

10  ps.  Hence,  plasma  displays  are  the  first  choice  at  this  time. 

Their  future  development  potential  is,  unfortunately,  not  good. 

2 

l:or  reasons  of  plasma  physics,  a  display  with  10000  points  would 

22" 

have  to  be  about  5  in  large.  The  high  driving  voltage  of  about 
1 2 S  V  is  also  undesirable.  Hence,  even  though  the  equipment 
liscussed  in  this  report  was  designed  for  a  plasma  display,  care 
vas  taken  that  most  of  the  electronic  circuits  could  be  used  with 
my  other  crossbar  display  too. 

The  crossbar  devices  with  the  host  potential  for  future  de¬ 
velopment  appear  to  he  those  based  on  acousto-opt ical  effects. 
There  are  a  number  of  reasons  for  this,  first,  the  area  of  a  point 
typically  a  light  transmitting  or  reflecting  spot  -  can  be  made 

m  th<*  order  of  10“  pm1"  as  compared  with  500“  pm  for  the  plasma 

2 

lisplays.  This  permits  potentially  displays  for  10000  points  to 
lave  manageable  size.  Second,  the  speed  of  such  devices  permits 
>nc  to  expect  a  transi/tion  time  from  one  pattern  to  the  next  in 
i  ps  or  less,  which  is  faster  than  the  plasma  display  can  achieve. 
One  should  note  that  such  acousto-opt i cal  scanners  exist  currently 
on  paper  only. 

from  the  standpoint  of  technology  one  should  thus  develop 
mage  analysis  with  the  help  of  the  available  plasma  displays, 

>ut  keep  in  mind  that  circuits  and  methods  should  be  applicable 
to  acousto-opt ical  crossbar  devices  too.  At  t ho  same  time,  image 
maly-is  provides  a  large  application  area  for  acousto-opt ical 
levices,  and  may  thus  help  stimulate  the  development  of  acousto- 
jptical  crossbar  scanners. 
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l'ig.  3.  Characterization  of  shape  by  u  listing  of  relative  li  ngths 
of  sides  and  angles,  if  two  sides  arc  equally  long  and  at  the  same 
time  longer  than  all  the  other  sides. 
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fig.  4.  Character i  .at  ion  of  the  shape  of  straight  line  structures 
t  iat  arc  not  closed. 
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l  it?*  5.  definition  of  curvature  (a),  and  two  possible  definitions 
of  a  radius  of  curvature  (b,c). 
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big  •  6  An  ellipse  on  a  field  with  a  finite  number  of  elemental 
points  or  squares  with  a  finite  area  Ax. \y  rathei  than  in  a  field 
with  nondonumera  hi  y  many  squares  with  infinitesimal  area  dxdy. 

The  crossed  squares  represent  the  ellipse,  the  solid  line  with 
the  no  re  familiar  shape  of  an  ellipse  is  only  shown  as  a  visual 
aide.  The  51  x  26  elemental  points  or  squares  represent  a  section 
of  a  liqtid  crystal  or  plasma  tube  display. 
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l;iK-  7.  T  bo  16  possible  contour  operators  produced  by  an  illum¬ 
inated  scanning  pattern  with  2x2  points. 
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Fig.  8.  Two-dimensional  Walsh  functions  wal(i,t,)  wal  and 

block  pulses  or  block  functions  blo(j,t;)  blo(k.T));  white  - 
illuminated  <  1  for  both  systems  of  functions;  black  »  -1  for 
the  Walsh  functions  but  black  =  opaque  *  0  for  the  block  functions. 
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F i  5 .  1.1.  The  20  sequences  of  three  points  complementing  the  20 
se jucnc  s  of  Fig.  10  by  numbering  the  points  in  the  opposite 
di  ecti  >n.  The  numbering  starts  with  the  point  denoted  0. 


As/Ax-1  |  As/A»-/2 

F i j*. •  12.  Tlic  throe  possible  values  of  the  length  element 
As(n)/Ax  in  point  n  on  a  square  field. 


Fig.  IS.  A  rectangle  on  a  square  field. 


Fig.  14.  lill  ipse  represented  by  a  discontinuous  set  of  a  finite 
number  of  points  (a)  and  by  a  continuous  set  of  nondenumorably 
many  points  with  a  discontinuous  first  derivative. 
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I- ]  g .  15.  LI  1  ipse  represented  by  straight  lines  between  coiner 
points.  The  angles  at  the  corner  points  are  always  non-negative 
which  means  either  positive  or  zero. 
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l-'ig.  16.  Typical  contour  of  an  opaque  area  (a),  definition  of 
the  angles  f(n)  and  At'(n)  of  the  straight  lines  between  corner 
points  (b) ,  and  t  lie  characterization  of  the  corner  points  by 
the  positive  steps  of  y(n)  (c). 
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Hig.  17.  Determin.-i  t  ion  of  r(n)  and  As(n)  from  the  values  X  and  Y 
of  a  cartesian  coordinate  system. 


)  ig.  18.  Rectangle  1  opresent ed  by  straight  lines  between  corner 
poii  ts.  1  he  angles  at  the  corner  points  are  always  non-negative 
hut  the  corner  points  5,  8,  11,  ...  are  degenerated  since  Ac 
equ:  Is  7. cro  there. 


Fig.  20.  Characteristic  values  of  the  ellipse  of  Pig.  6  as 
functions  of  the  normalized  curve  length  s*  from  the  point  n 


Fig.  21.  Total  rotation  of  the  tangent  around  a  closed  curve 
must  equal  2it ,  since  in  a  description  with  polar  coordinates  the 
variable  r  would  rotate  2rc  for  a  closed  curve,  and  the  tangent 
must  rotate  equal ly  much  to  make  the  tangents  at  the  beginning 
and  end  of  the  curve  coincide. 
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Fig.  22.  Use  of  the  normalized  distances 
along  the  straight  solid  lines  in  Fig.  IS 
corresponding  points  sJc*  and  s2c*  on  the 


Sj*  and  s2*  measured 
to  determine  the 
curved  line. 
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Tig.  23.  First  unsatisfactory  method  to  generalize  the  circular 
section  of  Fig.  19  to  eliminate  all  straight  lines  or  "flat  spots 


P(n*l) 


Fig.  24.  Second  unsatisfactory  method  to  generalize  Fig  19  i 
eliminate  all  fljt  spots. 


Fig.  25.  Satisfactory  generalization  of  Fig.  19  to  eliminate 
flat  spots. 


p(n.d)  *  p(n)  •  d  As"(n)/(  As(  n)  ♦  AsXn)l  I  a  n  Ae  (n) 
d ■  S-Sjln),  As(n)»As(n-1).  S[,(n).s(n)-As/n  1V2  for  As(n)«As(n-ll 
d*  s0(n)-s .  As(n)«As(n).  s^n)-s(n)*As(n)/2  lor  As(n-l)«As(n) 

As'(n).llAs(n)-As(n-1)|/2l 

lig.  26.  Definition  of  u*(s*)  for  any  value  of  s*  from 
X*(n)  defined  for  certain  values  of  s*  =  s*(n)  only.  (a.) 

Step  function  appi ox imav ion ;  (b)  Linear  interpolation  approxima¬ 
tion. 


81 


aaaaBaaaaaiaaaa  vaaaa  >>  aaaaaaaaaaa 

aaaaaaaaaaaaaaa>>*-';;.:.v.-..:  J»'.aai.'aaaaa 

aa9gaaiaaa.‘iic<-::5*3**&>!oc<-:vi,***a*M 

BBBBBBaBBB^;.**>>;«*S*c<c«******BBBBa 
a  a  a  a  a  a  a  a  a  ■;>  *  *  •:<  *  •:<  •:<  *  *  •:«  *  •:<  *  •:<  *  *  a  a  a  a 

a  a  a  a  -  -  a  .•  •!«  *  ■!<  *  *  *  {•  '5-  *  a  *  *  *  *  .?<  A  &  A  A  a  a 

a  I  *  >5  tf  s  C*  5  v1  i"  **■  v  $  o  *  *  A  {« >}  A  5  A  •;«  A  A  A  A  A  v  a  a  a 

a  a  c  A  *  •!<  A  A  >:•  A  A  •!■  A  «;•  A  •;  A  A  A  A  >:<>:>:•  a  a  a  a  a  a  a 
a  a  a  aA  >:«:<  A  A  A  A  *  A  A  A  •;>  A  *  A  A  A  c<  :<  a  a  a  a  a  a  a  a  a 
■■■■■■.•■"J  v  >:•  A  A  o* a i a  a  a  a  a  a  a  a 
a  a  a  a  ■  .  >  a  *  A  *  A  A  A  A  A  A  >;-  v  •;<  a  a  a  a  a  a  a  a  a  a  a  a 

BBaBaaai*v*  A-?--:-£-:<*ivBaaaeaBBBaaaBB 
aBaaBaamiaaakiaaaBaaaaaaaaaBaaaaa 


The  ellipse  of  Fig.  C  lotatcd  through  the  angle 


Fig.  29.  He  el  ipse  of  lig.  28  represented  by  straight  lines 
between  corner  p  ini'.. 


alized  curve  length  s*  from  the  point 


mbs 

iggHSaB: 


mm* 


l'ig.  31.  Connection  between  the  angle  of  rotation  9(b)  between 
two  ellipses  and  the  distance  between  their  main  axes  measured 
alo  ig  either  one  of  the  ellipses. 
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A  circular  curve  with  76  points  on  a  square  t'ivld. 


Fig.  37.  The  contour  of  the  circular  area  of  Fig.  36  or  the 
circular  curve  of  Fig.  .36  represented  by  straight  liner,  betw 
corner  points. 
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Fig.  38.  Characteristic  values  of  the  circles  of  Figs.  3(>  and 
37  as  function  of  the  normalized  curve  length  s*  from  the  point 
n=l.  Note  that  the  linear  approximation  of  Fig.  26b  is  used 
foi  K*/2  rather  than  the  staircase  approximation  of  Fig.  26a. 


0  10  20  30  40  50  60  70  80  90  100 

s’- 


Fig.  39.  Normalized  curvature  of  the  circles  of  Figs.  36  and 
3  ,  averaged  over  11  samples,  less  the  normalized  curvature 
2u/i00  of  an  ideal  circle  with  100  points  on  its  contour. 


plasma  display 


Fig.  40.  Block  diagram  of  the  scanning  device  for  contour  recog 
nition  in  aerial  photographs.  DDC  display  drive  control,  BFG 
binary  function  generator;  FSW  format  switches;  SSW  selector 
switches;  PUR  pown  drivers. 


